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ABSTRACT 


The  general  problem  chosen  for  consideration  concerns  the 
establishment  of  a  payload  in  o:~bit  about  the  planet  Mars.     Empha- 
sis in  this  studv  has  been  placed  upon  a  comparison  of  the  various 
paths  available  for  a  free   fall     space  stabilized  rocket,    the  re- 
quirements necessary  to  achieve  such  paths,   and  the  selection 
of  the  optimum  path,     The  optimum  path  has  been  compared  with 
other  trajectories  and  conclusions  have  been  drawn  as  to  the 
requirements  for   deviations  from  the  optimum  path.     The  study 
is  primarily  concerned  with  the  midcourse  portion  of  the  trip, 
which  has  been  defined  by  the  authors  as  that  part  of  the  overall 
trajectory  extending  from  a  desired  circular  orbit  about  Earth 
to  a  desired  circular  orbit  about  Mars.     A  firm  basis  for  the 
se  ection  of  the  orbit  has  beer,  established      Treatment  has  also 
been  given  in  this  report  to  the  ro-  tinuous  determination  of 
present  position  in  space  and  *be  initiation  of  corrective  action 


when  deviation  from  a  desired  path  exists.     The  study  also 
includes  generalized  instrumentation  and  guidance  requirements. 
Extensive  quantitative  information  has  been  complied  which  will 
serve  as  a  basis  for  the  determination  of  instrumentation  re- 
quirements as  well  as  creata  a  proper  background  for  an 
appreciation  of  flight  within  the  solar  system. 

Thesis  Supervisor:     Winston  R.    Markey 

Title:  Assistant  Professor 

of  Aeronautical 
Engineering 
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OBJECT 


The  object  of  this  paper  is  to  investigate  the  mid-course 
guidance  requirements  to  place  a  payload  in  orbit  about  Mars 
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CHAPTER  I 
INTRODUCTION 


The  authors  have  undertaken  this  study  in  order  to  bring 
together,    in  proper  context,    some  considerations  which  may  be 
anticipated  in  interplanetary  flight .     A  journey  to  the  planet  Mars 
has  been  chosen  as  the  basis  for  the  presentation.     However, 
many  of  the  ideas  and  methods  set  forth  here  can  be  considered 
applicable  to  flight  to  other  bodies  within  the  solar  system.     It 
is  felt  that  no  loss  of  generality  will  result  if  attention  is  focused 
upon  this  particular  choice. 

This  paper  has  evolved  as  a  study  of  several  problems  as- 
sociated with  interplanetary  flight.     The  objective  to  be  pursued 
is  concerned  with  placing  a  payload  in  orbit  about  the  planet  Mars . 
Departure  is  assumed  to  originate  from  orbit  about  Earth.     The 
attainment  or  possibility  of  flight  is  not  questioned.     Rather, 
this  is  a  broad  study  of  some  preferred  methods,    minimum  re- 
quirements,   and  ideas  concerned  with  fulfillment  of  the  general 
objective , 

Recent  developments  have  justified  studies  of  interplanetary 
flight.     It  is  realized,    however     that  the  topic  is  hopefully  opt:.- 
mistic  in  many  respects.     For  this  reason  the  study  is  neces- 
sarily broad  and  is  not  designed  to  give  an  intricate,    detailed 
account  of  methodology 

A  flight  to  the  moon  will  logically  precede  planetary  flight, 
but  its  close  proximity     causes  this  problem  to  become  quite 
specialized  as  contrasted  to  the  Mars  flight.     Piloting  a  boat  in 
an  inland  waterway  as  compared  to  navigating  a  ship  on  the  open 
sea  is  an  appropriate  analogy      The  two  problems  are  independent 
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in  a  large  sense. 

This  study  has  been  preceded  by  an  extensive  search  through 
the  literature  for  information  regarding  this  comparatively  new 
field.     Treatments  have  generally  varied  from  the  popularized 
texts  to  the  highly  specialized  mathematical  developments. 
There  appears  to  be  a  paucity  of  practical  approaches  to  the 
problem  which  can  be  interpreted  by  the  systems  designer  or 
engineer. 

The  authors  generally  envisage  an  unmanned  vehicle  with 
small  payload,   unpowered  except  for  terminal  and  corrective 
thrust  requirements .     It  is  believed  instrumentation  can  be 
handled  with  components  whose  performances  do  not  differ  too 
radically  from  slight  extensions  of  present  day  development. 

A  practical  simplified  view  has  been  taken  as  a  rule  due 
to  the  lack  of  an  available  high  speed  digital  computer.     Such 
a  device  would  of  necessity  be  required  if  accurate  trajectory 
studies  were  in  order.     However,   it  is  felt  that  little  departure 
from  reality  exists  under  the  assumptions  imposed.     In  fact, 
the  simplifications  add  to  the  comprehension  of  the  material. 

Solutions  are  sought  in  answer  to  the  following  questions: 

(1)  What  are  the  major  influences  of  the  solar  system 
upon  the  rocket  in  flight? 

(2)  What  simplifying  assumptions  can  be  made? 

(3)  What  routes  may  be  taken? 

(4)  What  criteria  are  available  for  the  selection  of 
one  trajectory  or  route  in  preference  to  another? 

(5)  How  may  position  in  space  be  continuously 
determined  for  navigational  purposes? 

(6)  How  may  the  rocket  be  assured  of  arriving  in  the 
target  area  if  errors  do  occur? 
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CHAPTER  2 
THE  SOLAR  SYSTEM 


A  summary  of  currently  accepted  data  concerning  the 
universe,    which  is  pertinant  to  this  study,    is  needed  for  formu- 
lation of  assumptions  and  as  basis  for  computation.     This  is 

presented  below  as  discussion  and  in  tabular  form  and  figures, 

2    4   9   44   53    54    59* 
and  has  been  derived  from  references.  '    '    ' 


2. 1       The  Sun  and  its  Planets 

The  solar  system  consists  of  the  sun,    nine  conspicuous 

planets,    and  numerous  smaller  planets.     The  conspicuous  planets 

are  placed  in  two  groups.     The  major  planets,    Jupiter,   Saturn, 

Uranus,    and  Neptune,    are  large,   have  low  density,    and  are 

surrounded  by  high  atmospheres.     The  terrestrial  or  minor 

planets,    Mercury,   Venus,    Earth,    Mars,    and  Pluto  are  similar 

to  the  earth  in  many  respects.     The  so  called  navigational  planets 

are  Venus,    Mars,    Jupiter,    and  Saturn,    deriving  this  classi- 

9 
fication  from  their  nautical  use  for  position  determination. 

For  the  period  of  1956-1970,   the  angular  motion  of  the 
major  planets  in  the  solar  system  will  be  as  shown  in  Fig.  2.1    . 
The  minor  planets,    Mars,    Earth,    and  Venus,   are  shown  in  their 
positions  as  of  1  January  1970,   and  in  their  respective  orbits, 
drawn  as  circles. 

A  generalization  which  serves  to  illustrate  the  area  occupied 
by  the  planets  and  sun  is  to  picture  a  piece  of  1/4  inch  plywood  cut 
in  a  circle  of  5  feet  diameter.     This  would  contain  all  the  planets 
in  their  orbits  around  the  sun  (except  Pluto)  if  drawn  to  a  scale 
such  that  Neptune's  major  planetary  axis  was  5  feet. 
*  Superscript  numbers  refer  to  bibliography. 
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Fig.  2.1       Angular  motion  of  major  planets  from  1965  to  1970 
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2.  2       Solar  Radiation  and  Temperature  in  Space 

The  solar  constant  of  sun's  radiation  is  1.  938  calories  per 

54 
minute  per  square  centimeter,  ±.05  This  value  is  equiva- 

lent  to  1.  35  x  10     ergs  per  second  per  centimeter  square  or 

87 . 1  x  10     ergs  per  second  per  inch  square  at  a  distance  of  1  AU 

from  the  sun.     This  incident  energy  varies  as  the  reciprocal  of 

the  range  squared,    as  shown  mFig.    2.  2 

The  maximum  temperature  of  the  surface  of  a  black  body 
planet  with  no  atmosphere  varies  as  the  reciprocal  of  the  square 
root  of  the  distance  of  the  body  from  the  sun,    as  shown  in  fig.  2.  3- 
At    a  distance  of  1  AU  from  the  sun  this  temperature  is  392   K. 


2,  3      Approximate  Orbital  Paths  of  Planets 

The  path  of  a  planet  or  body  in  orbit  about  the  sun  can  be 
completely  described  by  seven  elements,  in  so  far  as  the  body 
depends  only  upon  the  sun's  attractive  force.     These  are: 

1.  Semi  major  axis  (a) 

2  .  Eccentricity  (e) 

3.  Inclination  (i) 

4.  Longitude  of  ascending  node  (-^V) 

5.  Longitude  of  perihelion  (  n) 

6.  Period  (T) 
7  Epoch 

These  are  shown  in  Fig.  2.  4  and  2.5. 

Kepler's  Three  Laws  of  Planetary  Motion  are  also  appli- 
cable and  are: 

1.  The  planets  move  in  ellipses,   with  sun  at  one  of 
the  focuses 

2 ,  The  radius  vector  sweeps  out  equal  areas  in  equal 

time 
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Fig.    2.4      Basic  ellipse 
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Fig.    2.5     Basic  ellipse  oriented  to  reference  plane 
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3.   The  square  of  the  time  of  revolution  is  proportional 
to  the  cube  of  the  mean  distance. 

2.4      Perturbations  of  Planetary  Orbits 

The  actual  motion  of  the  planets  may  best  be  described  as 
follows: 

1.  The  orbits  of  the  planets  are  not  fixed  but  gradually 
change  in  eccentricity,    inclination,  (etc. ).     These 
changes  are  classified  as  secular  perturbations; 
however,   these  perturbations  are  so  slow  that  during 
any  one  revolution  the  orbit  is  altered  very  little  and 
during  the  time  of  any  given  single  interplanetary 
flight,    can  be  neglected. 

2 .  A  planet  in  its  motion  does  not  exactly  follow 
Kepler's  Laws  but  is  seen  to  vary  in  three  dimensions 
relative  to  its  ideal  Keplerian  Ellipse.     Such  changes 
are  very  small  and  are  called  periodic  perturbation. 


These  values  expressed  in 

heliocentric  longitude  (X  ) 

are: 

Mercury 

15" 

Venus 

30" 

Earth 

1' 

Mars 

2' 

Jupiter 

30' 

Saturn 

70' 

Uranus 

60' 

Neptune 

35' 

In  general  the  larger  the  terms 

the  longer  the  period  of 

periodic  perturbation.     For  example 

Jupiter  and  Saturn  exhibit 

periods  of  periodic  perturbations  of  about  913  years 
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2.  5      The  Invariable  Plane 

The  most  stable  property  of  the  solar  system  is  the  so 
called  invariable  plane.     This  is  defined  by  the  condition  that  the 
total  angular  momentum  of  the  solar  system  is  zero  about  any 
axis  in  this  plane.     The  invariable  plane  can  be  considered  as 
the  fundamental  reference  plane  of  the  solar  system  and  is 
inclined  1     38'  7"  to  the  ecliptic  reference  plane  of  1900  with 
ascending  node  in  longitude  106     42'  measured  fromT*   .     This 
inclination  is  decreasing  at  the  rate  of  .144"  per  year  and  will 
reach  47'  in  about  20,  000  years  at  which  time  it  will  begin  to 
increase  again. 

2.  6       Cislunar  Motion 

Cislunar  (earth-moon)  motion  in  the  plane  of  the  ecliptic  is 
that  of  two  unequal  revolving  masses  in  a  central  force  field, 
whose  center  of  mass  follows  a  slightly  distorted  Keplerian 
ellipse  of  eccentricity  (e)=.  01674.    This  cislunar  center  of  mass  is 
approximately  2880  miles  from  the  earth  center  in  the  direction 
of  the  moon.     See  Fig.  2.  6  and  2.  7, 

The  moon's  path  around  the  earth  has  an  eccentricity  .055 
and  a  mean  inclination  to  the  ecliptic  of  5     9' .     The  average 
distance  from  earth  to  moon  is  238,.  857  miles  (60.207  earth  radii) 
with  maximum  distance  of  252,  710  miles  and  minimum  distance  of 
221.,  463  miles. 

If  we  represent  the  cislunar  orbit  by  a  circle  of  radius  100 
feet,  the  moon  will  deviate  from  this  circle  approximately  3  inches 
on  either  side  in  one  transit  around  the  earth.     During  this  time 
the  geocenter  will  deviate  approximately  .04  inches  either  side 
of  the  orbit. 

2.  7       Solar  Motion 

Campbell  and  Moore  in  1920  found  the  solar  system  to  be 
moving  at  about  12.2  miles  per  second  with  respect  to  the  nearer 
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few  hundred  stars  toward  the  constellation  Hercules  ,     Jay, 
Sanford,    Stromberg,    and  Hubble,    among  others,    estimate  that 
the  solar  system  and  the  neighboring  stars  are  moving  at  about 
170  miles  per  second  toward  the  constellation  Cygnus . 

Conclusions  reached  by  Oort  place  the  solar  system  about 
32,  000  light  years  from  the  center  of  the  galaxy  which  is  in 
Sagittarius  and  that  this  position  is  about  two  thirds  of  the  dis- 
tance from  center  to  circumference.     At  170  miles  per  second 
the  period  of  galactic  revolution  of  our  sun  is  approximately 
2,  200,  000  years,   indicating  that  during  the  time  covered  by  any 
given  flight  within  the  solar  system,   the  sun  may  be  considered 
as  a  rather  stable  reference  point. 

2.8      Fixed  Stars 

The  nearest  star,  Alpha  Centauri,    is  at  a  distance  from  our 

13 
sun  of  4.3  light  years  or  2.  53  x  10      miles.     It  has  also  been  de- 
termined that  no  star  exhibits  a  parallax  of  more  than  .  76  seconds 
and  only  200  stars  have  proper  motion  greater  than  one  tenth  of 
"Barnard's  Runaway  Star"  which  traverses  10.25  seconds  of  arc 
per  year.     For  the  purposes  of  this  paper  a  "Fixed"  star,   as  the 
term  is  used  for  establishment  of  celestial  inertial  reference, 
will  be  considered  to  be  a  star  exhibiting  a  parallax  and  a  proper 
motion  of  less  than  one  quarter  second  of  arc.  Of  the  stars  checked 
to  date  only  about  .  01%  do  not  fall  within  this  classification  of 
"Fixed"  stars  . 

2.  9     Gravitational  Effects 

Newton's  general  Law  of  Gravitation  states  that, "Any  two 

particles  of  matter  attract  each  other  with  a  force  proportional 

to  the  product  of  their  masses  and  inversly  proportional  to  the 

square  of  the  distance  between  them"       or  f  =  K — 5 — ,   where  K  is 

BT 
the  so  called  Gaussian  gravitational  constant  equal  to 
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-8             cm3 
6.67  x  10"      s"    •     From  this  the  acceleration  of  gravity 

gM(m)  sec 
is  then  g  =  ^M.  . 

IT 

Values  of  the  acceleration  of  gravity  for  the  various  solar 

bodies  are  tabulated  in  Table  2.  2      with  the  mean  surface 

gravity  for  earth  equal  to  1.      Fig.    2.8    presents  "g"  versus  the 

radial  distance  from  the  earth,   mon  and  Mars  and  the  variation 

in  the  effect  of  the  sun's  gravitational  attraction  over  the  distance 

between  the  orbits  of  Earth  and  Mars . 

18 
Herrick  has  pointed  out      that  Gauss  was  able  to  determine 

2 

mathematically  the  value  of  K  ,   for  heliocentric  orbits  with  dis- 
tance in  AU,   accurate  to  nine  significant  figures  but  neither  the 
laboratory  value  of  gravitational  acceleration  nor  laboratory  units 
of  measurements  at  present  will  permit  any  greater  accuracy  than 
3  or  4  significant  figures. 

An  interesting  discussion  of  Gravitational  and  tidal  effects 

which  has  direct  bearing  on  any  travel  between  Earth  and  Mars 

56 
has  been  presented  by  Schaub      .      Schaub  states  that  the  inner 

54 
most  moon  of  Mars,   Phobos,    is  within  Rocke's  limit      .     This 

means  that  since  Phobos  is  assumed  to  have  the  same  density  as 

Mars  and  is  traveling  in  an  orbit  of  less  than  2.  44  times  Mars 

radius  (see  Table  2.2  )  its  surface  is  now  free  of  all  loose 

material.     Eventually,    due  to  the  bombardment  by  rays  from 

space,    more  and  more  of  this  moon  will  be  torn  away  forming  a 

ring  of  dust  particles  around  Mars .     This  implies  that  any 

terminal  orbit  about  Mars  should  stay  well  clear  of  Phobos  less  it 

hasten  this  disintegration  possibly  destroying  the  vehicle  in  the 

process. 
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CHAPTER  3 
THE  PROBLEM  SIMPLIFIED 

3.1       Introduction 

The  Law  of  Universal  Gravitation  as  enunciated  by  Newton 

states  "Every  particle  of  matter  in  the  universe  attracts  every 

other  particle  with  a  force  which  acts  in  a  line  joining  them,    and 

whose  intensity  varies  as  the  product  of  the  masses  and  inversely 

as  the  squares  of  their  distances  apart." 

44 
Moulton       presents  an  excellent  discussion  of  this  law  and 

states,  "it  will  be  observed  that  the  law  of  gravitation  involves 

considerably  more  than  can  be  derived  from  Kepler's  Laws  of 

Motion;  and  it  was  by  a  master  stroke  of  genius  that  Newton 

grasped  it  in  its  immense  generality  and  stated  it  so  exactly  that 

it  has  not  been  necessary  to  change  a  single  syllable  in  over  two 

hundred  years.     When  contemplated  in  its  entirety  it  is  one  of  the 

grandest    conceptions  in  the  physical  sciences.  " 

It  is  evident  that  the  various  bodies  comprising  the  Solar 
System  will  exert  certain  forces  upon  a  rocket  in  flight  among 
them.     Because  of  the  extremely  large  masses  and  distances 
involved  we  will  omit  consideration  of  the  minor  bodies,    such  as 
asteroids  and  meteor  swarms,   and  show  later  by  inference  that 
this  omission  is  justifiable „     Considering  then  only  the  Sun,    Moon, 
rocket,    and  nine  planets     we  are  faced  with  a  twelve -body  problem 
It  would  appear  expedient  to  simplify  this  problem  if  possible 
since  twelve  bodies  require  thirty -six  second  order  differential 
equations  (thirty -three,    if  referred  to  the  Sun)  to  adequately  de- 
scribe their  motions. 


3.  2      Differential  Equation  of  Motion  of  the  Rocket 

In  order  to  reduce  the  complexity  of  the  problem  let  us  first 
examine  the  terms  of  the  equation  of  motion  of  the  rocket  and 

consider  their  magnitudes.     With  reference  to, Fig.  3.1   define 

71 
I  as  an  inertial  reference  frame  in  the  convention  of  Wrigley 

and  0  as  a  frame  within  the  Solar  System  non-rotating  with  respect 

to  I  and  with  axes  parallel  to  those  of  I.     Note  that  at  this  point  0 

is  restricted  only  to  lie  somewhere  in  the  vicinity  of  the  Solar 

System.      A  represents  the  position  of  the  rocket,    a  variable. 

We  may  write  Newton's  Second  Law  in  vector  notation  as 
follows; 


t*Wl  "5^       ■     ff    +    fnf  <3-» 

where  [  RT  A]  T  is  the  second  derivative  of  the  vector  RTA  with 
respect  to  the  frame  I  and  the  other  terms  are  previously  defined, 


Then; 

Tf    =    ^©A    +  %A    +  ^(TA  +  5?A    +    ^OA  +  %A 

+  54A  +  ^hA  +  %A    +    5*A    +    GEA  (3-2) 

where  G.    .      is  defined  in  Derivation  Summary  1  . 

If  there  is  a  particle  at  point  0  of  mass  m  '  then; 

mJJ*TJT  =  m     (CU     +  Gm      +   G-      +   G  w ) 

ol     IoJ  I  o      0o  ©o  do  Qo 

(3-3) 
and  the  masses  may  be  cancelled. 

By  geometry 

["iJlMloJi  +  lSjj!  (3-4) 


Fig.    3.1      Reference  frames 


ORBIT  OF 
VENUS 


ORBIT  OF 
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(a)  INFERIOR   BODY 


(b)  SUPERIOR  BODY 


Fig.    3.2      Planetary  configurations  for  Table  3.1 
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Hence  by  substitution 

t%oll  +  lK0A]l=Tf+Tnf  (3"5) 


or 


(G©o+G®o+G<ro -,+  tRoA]rfnf 

+<G©A  +  G®A  +  G«A  + >  <3-6) 

Hence  in  full 

[  ^o  Al  I  ■  Tnf  +  <5©A  -  W  +  (G®A  -  "W  +  (GCA  "  Gjo> 

+(g?a  -%)  +  (V "  V+(e<& " 5  °"o) 

+<G2)A  "  5*0>  +  <BhA  "  Gflo)+  (5©A-  %o> 
+<Sa"  V)  +  (GEA-Geo»  (3-7) 

We  shall  now  consider  applying  this  derivation  (frame  0)  to 
various  members  of  the  Solar  System. 

3,  3      Inertial  Points  Within  the  Solar  System 

Looking  first  at  the  motions  of  the  Sun,   or  the  Solar  System 
as  a  unit,   it  appears  from  a  statistical  analysis  of  the  motions 
of  the  "nearer"  stars  that  the  Solar  System  possesses  a  velocity 
of  approximately  12.2  miles  per  second  relative  to  these  "nearer" 
stars.     The  "nearer"  stars  range  from  4  to  300  light-years 
distant  from  the  Sun.     Further,   the  Sun  is  thought  to  occupy  a 
position  approximately  20,  000  light -years  from  the  center  of  our 
Galaxy  (The  Milky  Way)  with  a  period  of  revolution  about  the 
galactic  center  on  the  order  of  2,  200,  000  years.      (See  Section  2.  7 
for  further  details  on  the  Solar  System).     In  addition  our  Galaxy 
is  a  member  of  the  so-called  Local  Group  consisting  of  17  galax- 
ies,   clusters,    and  clouds  including  the  Andromeda  Galaxy,   The 
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Large  and  Small  Magellanic  Clouds,    and  several  small  clusters. 
It  is  thought  that  these  17  members  rotate  about  a  common  center 
of  mass  located  somewhere  between  our  galaxy  and  the  Andromeda 
Galaxy.     Lastly,    current  theories  on  the  "Red  Shift"  of  the  optical 
spectra  of  distant  galaxies  conclude  that  these  galaxies  are  re- 
ceding from  the  Solar  System  at  such  tremendous  relative  ve- 
locities that  certain  effects  attributed  to  the  Theory  of  Relativity 
must  surely  be  present.     It  is  not  the  purpose  of  this  paper  to 
argue  any  of  the  current  theories  or  to  attempt  any  analysis  in- 
cluding Relativity  Effects.     These  motions  to  which  the  Solar 
System  is  subject  are  mentioned  solely  as  a  basis  for  the  following 
statements: 

The  Sun  is  probably  undergoing  motions  involving  intricate  and 
complex  accelerations;  however,   the  vast  macrocosmic  scale  of 
distances  involved  is  such  that  during  the  time  of  flight  of  our 
rocket  the  effect  of  these  ill-defined  accelerations  is  neither  ap- 
preciable nor  measurable.     This  then  permits  the  statement, 

"Newtonian  space  is  that  unaccelerated  with  respect  to  the  'fixed' 

71 
stars",       particularly  since  it  is  not  known  that  any  stars  are 

fixed,   but  only  that  their  motions  are  very  small  during  recorded 

history. 

The  foregoing  definition  of  Newtonian  space  is  preferred  by 
the  authors  since  it  appears  at  least  semantically  more  defi- 
nitive than  the  statement  "Newtonian  space  is  that  space  wherein 

71 
Newton's  Laws  are  valid."        We  say  then  "The  possible  motions 

of  the  center  of  the  Solar  System,   however  complex,   within  our 

detective  ability  are  not  known  to  interfere  measurably  with  the 

information  obtained  while  ignoring  these  motions  and  effects 

thereof.  " 

3.  4      Simplified  Differential  Equations  of  Rocket  Motion 

On  the  foregoing  basis  we  allow  ourselves  to  define  the  center 
of  (1)  the  Sun,    (2)  the  Earth,   and  (3)  Mars  each  in  turn  as  the  lo- 
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cation  of  our  inertial  reference  frame  0,     Equation  (3-7)     may 
then  be  written  successively; 


fR©A]I  =  fnf  +  SOA+(<W%<£f(6<lA  ^«©  }  '  • 
f  ^®A]  I  =  7nf  +  S©A  +(SOA  "  %©  )+  {^(LA~  S(L® ) 
l*&h  =  ~nf  +^dA  +  (50A  "  GOCf>  +  (G(tA-  %j) 


(3-8) 
(3-9) 
(3-10) 


2 
Table   3.1     gives  values  in  ft/sec     of  the  maximum  values 

the  above  terms  in  parenthesis  will  hold,   using  mean  values  for 
the  orbital  radii  of  the  bodies.     Use  of  these  mean  radii  is  con- 
sidered justifiable  since  only  rarely  will  the  bodies  be  aligned 
such  that  these  maximum  values  actually  occur.     Figure   3.  2 
clarifies  the  manner  in  which  these  values  do  occur. 

By  inspection  of  Table  3.1;  if  we  establish  that  all  acceler- 
ations  less  than  10        ft /sec     are  to  be  ignored,   we  may  eliminate 
from  future  consideration  most  of  the  parenthetic  terms  in  the 
foregoing  equations.     We  then  have; 

•  • 

[ ^  ©A]  I  =  7nf  +  ^OA  +(^©A  ~  %D)  +  (^«A  "  <W 

+  (5d'A"5d,0)  (3"U) 

[lWl  =  7nf  +  %A  +(^OA  ^0©)  +  <G(IA"^(I©  ] 

+(5dA  "  W  (3"12) 

[*tfAl  I  =  Fnf  +  %A  +  (5OA-  ^x5>  +  <G«  a"  %cf) 

+(G©A"  <W)  (3-13) 

We  have  now  reduced  the  problem  to  one  of  five  bodies  rather 
than  twelve.     We  will  investigate  the  possibility  of  reducing  it  even 
further. 


t  -  2C_.  47  5  days 
4      2 


lloon 
Orbit 


1^ 


t  =  13.65  days 


R 


1-1 


240,  000  miles 


R2-2     =  563j  00° 
R33|  =  653,000 


|R4_4|  =  753,000        " 
Fig.  3.  3     Illustration  of  Rocket-Moon  Separation 


3.  5      The  Lunar  Gravity  Field  Effect 


22   . 


Referring  to  the  Hohmann  Ellipse   '    it  may  be  shown  that  in 
the  first  15.2  days  of  rocket  flight  along  this  trajectory  the  radius 
vector  of  the  rocket  measured  from  the  center  of  the   Earth 
increases  by  approximately  600   000  miles.     During  this  period 
the  rocket  is  in  almost  linear  flight,   hence  the  magnitude  of  the 
Lunar  gravity  field  is  easily  estimated  and  is  shown  in  Fig.  3.3. 

At  the  end  of  the  15.  2  days  the  value  of  the  gravity  field  has  de  = 

-5  2 

creased  to  less  than  1  x  10       ft/sec     and  is  continually  decreasing 

As  a  result  of  these  rudimentary  estimations  it  is  evident 
that  the  effect  of  the  lunar  gravity  field  need  be  considered  only 
during  the  initial  portion  of  the  trajectory.     The  example  selected 
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above  is  unique  in  that  it  has  been  shown  to  require  the  least  possi- 
ble   rocket  velocity  hence  the  rocket  would  be  subjected  to  Lunar 
influence  for  a  longer  time  than  in  any  other  trajectory.     The  fore- 
going estimates  then  are  an  upper  bound  to  the  values  which  actu- 
ally may  occur,  * 

3.  6      The  Flight-Phase  Concept 

Similar  arguments  to  those  of  the  preceding     section  may  be 
employed  to  justify  ignoring  the  effect  of  the  Earth  and  Mars  grav- 
ity fields  over  certain  portions  of  the  flight.     Certainly  if  a  com- 
puter is  designed  to  employ  the  values  of  these  fields  there  will 
be  a  certain  sensitivity  level  below  which  no  cognizance  of  these 
values  will  be  taken. 

Thus  we  may  write  with  some  degree  of  justification  the 
following  equations; 

[*©Jl=Inf    +    G©A  Phase  II  (3-14) 

[*Wl  =  Fnf    +   ^"©A  PhaSeI  <3-15> 

f%A]I"_fnf    +    GdA  Phase  III  (3-16) 

where  the  limits  of  the  Phases  are  to  be  identified. 

Thus  by  consideration  of  the  trajectory  and  the  magnitudes 
of  the  forces  involved  we  reduce  our  problem  from  one  of  12  bodies 
to  three  separate  two -body  problems. 


3d 


TABLE    3.1 


This  table  displays  values  for  the  parenthetic  terms  of 
equation  .3-  7     In  the  general  form  we  have; 


(G.A-G.   )l=  k.  R   .2(-^-t>    -    -i-s 
jA       jo  I        j      03  2  - 

jA  jo 


) 


where  k.  is  the  surface  gravity  of  the  indicated  body 

R    .  is  the  radius  of  the  body  in  AU 
oj  J 

R..  is  the  distance  separating  the  body  and  the  rocket  in  AU 
R.     is  the  distance  separating  the  body  and  point  0  in  AU 

J 


*T 


erm 


G0A~GO 

G©A"G© 
GdA~GG 


G 


G> 


gqa~go 

cfA-GCf 


G 


G2|A"G4 
GhA"Gh 

g6a"g6 

G  *A~G¥ 
G  EA-Gg_ 


0=0 


-32.2 

~4 

1x10  ^ 

2.025  x  10' 
1.635  x  10 
-12.25 
7.1  x  10" 
1.732  x  10' 


4.65  x  10 


-10 


1.585  x  10 


t6.62  x 10 


-10 


-13 


0  =  © 


2.78  x  10' 


*  1  x  10 
8.03  x  10' 
6.21  x  10" 

-12.25 
9.21  x  10" 


4.025  x  10' 


7.17  x  10 


2.01 x 10 
tl.01  x  10" 


■10 


■10 


•12 


o  =  cT 


2.78  x 10' 

-32.2 

*  1  x  10"' 

8.03  x  10' 

6.21  x  10" 


9.21  x 10 


-7 


4.025  x 10 
.-10 


8 


7.17  x  10 


2.01 x  10 


■10 


+1.01  x  10 


-12 


*  by  judicious  choice  of  firing  time 

+  Assuming  that  Pluto  and  Earth  have  the  same 

value  for  surface  gravity. 
+  Sample  calculations  are  given  in  Table  3.2. 


TABLE  3    2 
Sample  Calculation  for  the  Values  Listed  in  Table  3   1 
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CHAPTER  4 

THE  TRAJECTORY 
A.   SOME  POSSIBLE  TRAJECTORIES 

4 .,  1       Basic  Criterion  for  Selection 

Selection  of  the  type  of  trajectory  to  be  followed  during  the 
mid -course  phase  of  an  interplanetary  exploratory  trip  of  any 
type  is  dependent  upon  the  basic  criterion  of  when,   in  the  future 
we  contemplate  attempting  the  flight.     If  we  desire  to  make  the 
attempt  within  less  than  ten  years  it  is  reasonable  to  postulate 
that  presently  developed  power  plants  and  guidance  systems  must 
be  used  with  considerable  refinement  a  necessity.     If  the  attempt 
will  not  be  made  until  some  time  after  1970,   more  exotic  power 
plants  and  guidance  can  probably  be  postulated.     However,   the 
basic  trajectories  can  be  summarized  and  a  selection  made  for 
further  study  in  this  paper. 

4,  2     Classifications  of  Trajectories 

All  trajectories,   based  on  thrust  available,    can  be  divided 
into  three  classifications:  (1)  Powered,    (2)  Unpowered  and  (3) 
Combination  of  Powered  and  Unpowered.     Although  this  statement 
is  simple,   the  actual  engineering  and  instrumentation    re- 
quirements are  exceedingly  complex  and  somewhat     dissimilar. 

4.  3     Powered  Trajectories 

The  term  "Powered  Trajectory"  implies  continuous  thrust  with 
either  acceleration  or  deceleration  used  during  the  entire  flight. 
Also  implied  is  the  ability  to  apply  this  thrust  in  any  direction  by 
multiple  thrust  units  or  vehicle  reorientation  with  respect  to  the 
trajectory. 
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The  upper  limit  of  this  type  trajectory  can  be  considered  to 
be  strictly  interceptor,   while  the  lower  limit  is  that  of  a  distorted 

Keplerian  Ellipse .     Interceptor  trajectories,   using  proportion- 

43 
al  navigation,    collision  course,    line  of  sight,    or  constant  lead 

similar  to  present  day  antiaircraft  and  air-to-e.ir  guided  missiles, 
require  nearly  unlimited  thrust.     However,    many  simplifications 
in  guidance  methods  can  be  realized  due  to  reduced  time  of  flight 
and  reduction  in  accuracy  requirements  if  high  speed  inter- 
ceptor methods  are  used. 

With  small  thrust  power  units  such  as  the  ion  gun,    solar 

furnace,   or     Nuclear -chemical,    etc.  the  trajectory  can  be 

14 

a  logarithmic  spiral  such  as  developed  by  Forbes      or  a  series 

of  constantly  changing  ellipses  with  the  actual  trajectory  being 
the  envelope  oi  the  series.     In  addition  the  vehicle  could  be 
forced  to  traverse  a  predetermined  elliptical  path.     This  latter 
method  would  be  extremely  wasteful  of  both  fuel  and  time  , 

Powered  trajectories  are  difficult  to  compute  and  almost 
infinite  in  number.     The  practicality  of  detailed  powered 
trajectory  study  with  any  more  than  a  general  overall  look 
at  the  field  appears  to  be  limited  at  present  until  the  perfection  of 
suitable  thrust  devices  warrant  the  effort. 

Unpowered  Trajectories 

Basically  the  unpowered  trajectory  is  a  ballistic  flight  - 

22    43    72 
incorporating  the  Fire  Control  Problem      '     s'        with  that  of 

21   44   60 
determining  the  path  of  a  comet  in  space.      '       '         This  implies 

the  application  of  a  large  impulsive  thrust  applied  over  a  very 

short  interval  of  time  with  respect  to  total  time  of  flight.     This 

would  result  in  a  hyperbolic  path  with  respect  to  a  planet  during 

application  of  thrust.  J    '        At  the  end  of  this  time  the  hyperbolic 

trajectory  would  merge  with  an  elliptical  path  with  respect  to  the 

sun,   and  the  vehicle  would  be  in  free  flight,    acted  upon  basically 

by  gravitational  attraction  of  the  sun  with  slight  perturbations 
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caused  by  the  planets,    iseeSection  3.6), 

A  trajectory  of  this  sort  leads  to  simplification  or 
elimination  of  a  guidance  system  for  mid-course.  However., 
extreme  accuracy  in  the  orientation  of  the  velocity  vector  at 
termination  of  impulse  is  mandatory.  Very  small  errors  in 
magnitude  or  direction,  or  unforseen  gravitational  gradients 
during  transit  can  cause  large  errors  at  the  trajectory  terminus 
as  shown  in  Section  8.6. 

4„  5       Combination  Trajectories 

Trajectories  containing  thrust  periods  and  ballistic  no- 
thrust  periods  are  complicated  by  virtue  of  attempting  to  combine 
two  separately  difficult  problems  into  one.     This  method  does 
permit  slightly  lower  accuracies  of  the  initial  velocity  vector 
but  requires  accurate  mid-course  position  determination  to  be 
of  any  use. 

The  basic  advantage  of  this  method  is  that  presently  developed 
thrust  units  can  be  used  for  correcting  or  changing  the  trajectory 
during  mid-course,   thereby  increasing  the  probability  that  the 
mission  will  be  completed.     Using  contemplated  low  thrust  cor- 
rective  units      '        powered  phases  could  be  of  long  duration  and 
be  equally  as  effective  as  high  thrust,    short  burning  time  chemical 
units . 

4.6      Trajectory  Selected  for  Further  Study 

The  basic  mid-course  trajectory  chosen  for  further  study  in 
this  paper  is  that  of  a  combination  of  unpowered  and  powered. 
The  ideal  or  hoped  for  situation  is  one  in  which  no  corrective 
thrust  is  needed  during  the  time  of  transit.     This  would  mean  that 
a  relatively  unperturbed  ellipse  could  be  flown  within  the  limits  of 
terminal  accuracy  requirements.     However,   if  it  becomes  ap- 
parent during  transit  that  corrective  action  is  needed,   a  new 
ellipse  is  computed  and  proper  thrust  applied  to  get  on  the  new 
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trajectory.     The  position  determination  accuracy  requirements 
for  such  a  system  are  rather  high  as  are  the  thrust  control 
requirements.     However,    it  is  felt  that  these  requirements  are 
not  entirely  beyond  present  day  engineering  capabilities. 
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B,   THE  UNPOWERED  TRAJECTORY 


4.  7     Introduction 

To  intelligently  decide  upon  a  trajectory  to  be  traveled  to  Mars 
one  must  weigh  many  aspects  of  the  problem.     In  previous  por- 
tions of  this  paper  we  have  investigated  some  properties  which 
influence  a  vehicle  in  flight  in  the  Solar  System;  and.   in  addition, 
investigated  some  vehicular  behaviors  and  parameters  which 
could  be  selected  or  adjusted  to  our  benefit. 

As  previously  described,   the  "unpowered"  trajectory  is  a 
trajectory  provided  with  an  initial  impulsive  velocity,   which,   when 
perfectly  instituted,   needs  no  corrective  action.     This  type  tra- 
jectory lends  itself  well  to  analysis  since  the  vehicle  is  essentially 
a  body  acting  in  a  central  force  field  (subject  to  perturbations). 
Primarily  due  to  considerations  of  available  fuel  and  performance 
characteristics  it  was  decided  early  in  this  work  to  proceed  with 
an  analysis  of  this  type  trajectory.     It  is  emphasized  here  that 
this  type  trajectory  may  by  no  means  be  considered  more   at- 
tractive than  any  other  because  of  any  inherent  property  associ- 
ated with  it.     It  simply  appears  more  feasible  due  to  the  current 
state  of  the  art  in  the  development  of  rocket  vehicles. 

22 
Hohmann  "  appears  to  be  the  first  author  (circa  192  5)  to 

seriously  consider  the  solution  of  space  flight  problems.     One  of 
the  most  important  results  of  his  excellent  study  is  the  description 
of  the  minimum  energy  ellipse,   which  as  the  name  implies  re- 
quires the  least  rocket  velocity  of  any  configuration  of  two  planets. 
This  trajectory  has  become  known  as  the  "Hohmann"  ellipse  and 
has  the  characteristic  that  the  arrival  of  the  rocket  at  Mars 
occurs  when  Mars  is  diametrically  opposite  Earth  position  at 


time  of  firing  in  their  respective  orbits  about  the  Sun.     Note 
Fig.     4.1     . 
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Preston -Thomas     '  ^  J  agrees  with  Hohmann  and  shows  some 

simple  relations  of  the  Hohmann  ellipse  with  other  ellipses. 

Other  than  this  the  literature  shows  a  decided  lack  of  information 

on  other  trajectories  of  the  "unpowered"  type. 

It  is  the  purpose  of  this  chapter  to  show  some  relations  exist- 
ing among  the  "unpowered"  trajectories.     These  are  the  tra- 
jectories which  when  corrected  for  errors  are  described  as  the 
"combination"  trajectories  of  Section  4,5 


4.  8      Basic  Relations 

The  basic  relations  of  a  vehicle  in  motion  under  the  influence 

A4  60 

of  a  central  force  field  are  set  forth  by  Moulton      ,    Smart 

?  fi 
Kooy  and  Uytenbogaart      ,    and  others.     Those  general  equations 

selected  by  the  authors  which  appear  to  best  fit  the  problem  are: 

2, 


a(l-e   ) 

1-e  cos  6 


(4-1) 


V   =• 


(2/r  -  1/a) 


(4-2) 


,  e  sin  6 

tan  0    =  , a 

Y       1-e  cos  6 

47 
A  fourth  equation       which  will  be  of  interest  is 


(4-3) 


Jl 


e     sin©  ,     .    -1 

+  sin 


1-  e  cos  6 


e  -  cos  9 


1-e  cos  6 


+  2  T 


(4-4) 
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A  general  trajectory  situation  is  shown  in  plan  view  in  Fig. 4.  2 
Figs.    4.  2   and  4,  3  illustrate  the  quantities  appearing  in  the 
foregoing  equations  . 


Hohmann  Ellipse 


Fig.    4,1    The  Hohmann  Ellipse 


±.  9      Some  Simplified  Trajectory  Parameters 

We  logically  begin  our  analysis  by  considering  the  effect  of 
the  most  massive  member  of  the  Solar  System,   the  Sun,     Re- 
ferring to  Fig.  4    2  ;  for  a  given  trajectory  we  may  establish 
initial  and  finai  positional  conditions  (i.e,  -  we  know  we  are  on 
the  Earth  and  hopefully  we  will  arrive  at  Mars).     Hence  we  write 
(4-1)  as 


and 


r     = 
o 


a(l  -eT) 

1  -  e  cos  6 


o 


,      =     a(l  -  ez) 
A     1  -  e  cos  6 


A 


(4-5) 


(4-6) 


rocket  ellipse 
aphelion 


Fig.   4.2     General  t  rajectory  configuration 


Fig.   4.3      Some  vector-angle  relations  for  the  general 
trajectory 
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where  rQ   =   |  R0CJ  |  at    t  =  tQ,   and  rA  =  |  R0(5|at  t  -  tQ  +  tf . 
Similar  definitions  apply  to  the  angles     6     and  0.  . 

O  A 


Let  r     =   1  AU 
o 


Let  rA  =  1.5237  AU 


Dividing  (4-5)  by  (4-6)  and  rearranging,   we  have 


e 

0.5237 

1.5237  cos0A  - 

cos   6 
o 

Then  from  (4. 


(1  -  e  cos  6    )  r 

_, o      o 

I         2 
1  -  e 


(4-7) 


(4-8) 


which  may  be  expressed  more  conveniently 


1 

a 

1          2 

1  -  e 

1  -  e  cos  6 
0 

(4-9) 


where  a  is  in  AU. 

This  defines  two  of  the  parameters  (a  and  e)  of  any  ellipse. 

We  may  now  establish  the  vector  velocity  of  the  rocket  with 
respect  to  the  Sun  in  terms  of  the  trajectory  parameters. 


Le'Vo^     V€>A 


at    t  =  t    .     Then  (4-2)  may  be  written 


V     -  97.  500  *s/i  -  1/s 
o 


ft/sec 


where  a  is  in  AU 


Let  $    be  the  angle  shown  in  Fig.  4.4  at    t  =  t    .     Then 


(—                             — 1 

vo 

:   tan" 

e  sin  6 
0 

1  -  e  cos  6 
0 

' —                     — '      1 

(4-10) 


(4-11) 
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TO  SUN 


Fig.   4.4      Velocities  relative  to  Earth 


TO  SUN 


Fig.   4.5     Velocities  relative  to  Mars 
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Additionally,  from  Kepler's  Laws 

T-7m^  *3/2  <4=12) 

and  the  period  of  the  ellipse  is  known. 

The  boxed  equations  summarize  pertinent  parameters  in  the 
study  of  rocket  motion  with  respect  to  the  Sun,   and  form  a  basis 
for  the  calculations  of  velocity  vector  requirements  with  respect 
to  Earth  and  Mars . 

It  might  be  argued  here  that  greater  accuracy  should  have 
been  demanded  of  these  basic  equations  (i.e.    -  Use  1,523688  AU 

for    R  ^ft]     instead  of  1.  5237  AU).     The  astronomical  unit  is 

^  '  54  4 

known  accurately       only  to  about  one  part  in  10     hence  no  good 

purpose  could  be  served  by  extending  the  required  accuracy  of 
these  equations  beyond  the  fourth  decimal,    at  least  for  a  pre- 
liminary analysis . 

Calculations  were  carried  out  for  about  50  elliptical  tra- 
jectories and  the  values  for  e  and    Vqa    are  plotted  in  Fig     4.  6 
and    4.  7     as  aids  in  future  studies.     These  plots  presuppose  that 
6     and  0 A  are  given,   known,   or  selected.     A  logical  basis  for 

the  selection  of  6     and  6  .   is  the  subject  of  another  section, 
o  A  J 

Elhpticity  <  1  was  arbitrarily  chosen  since  excessively  high 
rocket  velocities  are  required  for  parabolic  or  hyperbolic  tra- 
jectories . 

4   10    Rocket  Velocity  Vectors  with  respect  to  Earth  and  Mars 

In  the  preceding     section  equations  were  developed  to  express 

parameters  of  the  elliptical  trajectory  and  initial  velocity  of  the 

rocket  with  respect  to  the  Sun.     Obviously,    since  Earth  and  Mars 

are  also  bodies  in  motion  with  respect  to  the  Sun  we  want  the 

velocity  of  the  rocket  with  respect  to  these  bodies .     The  relations 

are  simple,    geometric  in  nature,    and  are  developed  herewith. 

At    t   =    t 
o 
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v©A  =  vaA  +  voa 


(4-13) 


At  t  =  t     +  t- 
o         I 


V0A=VdA  +  V0d 


(4-14) 


as  defined  in  Definition  Summary    1 
From  equation  (4-3) 


6     =   tan 
Yo 


-1 


e  sin  6 


1  -  e  cos  6 


o 


And  from  Fig.     4,  4 


(4-15) 


V. 


®T: 


=  tan 


-1 


V     sin  6 
o ZSL 


V     cos  $      -  V/ 


o 


o 


'6 


(4-16) 


V, 


0A 


vo  cos  ^o  -  voe 

cos  \y0A 


The  properties  of  an  ellipse  identify 

r 


V 


-2-v 

rA       ° 


ifvA=  v0 


A       at     t  =  t     + 1- 
A  o        i 


<ji   .   =   tan 


-1 


e  sin  0 


A 


1  -  e  cos 


A 


and  from  Fig.  4.  5    and  equation  (4-14). 


%z 


=  tan 


VA  Sin  ^A 

Vl-'VAcos  ^A 


(4-17) 


(4-18) 


(4-19) 


(4-20) 


vcf< 


Vo6-  VA  cos  *A 


COS 


YdA 


(4=21) 


If  r     =  1.0000  AU  and  r.   =  1,5237  AU  then 
o  A 


V0©  =    97,  500  ft/sec 


and 


V0<5  =  78'  "°  ft/sec 
Figures  4.  8  and  4,9  show  V^  »  and  V^.  as  functions  of     6 
and  0».     These  values  are  important  primarily  as  a  step  in 
solving  the  trajectory  problem  and  are  given  here  for  possible 
future  use  by  the  reader. 

4.11    Time  of  Flight 

The  calculation  of  time  of  flight  is  necessary  at  this  point  in 
order  to  determine  the  angular  separation  of  Earth  and  Mars  at 
the  instant  of  rocket  firing.     Furthermore,   it  is  necessary  to 
know  U  in  sidereal  or  calendar  time  since  the  inclination  of  Mar's 
orbit  with  the  ecliptic  is  a  function  of  real  time. 

If  we  substitute  t     and  t,  as  defined  in  Definition  Summary  1 
each  in  turn  in  Equation  (4=4)      and  recognize  0A«»  0     and 

A  C 

6   EE  9,  we  have 
o        d 


jVl-e     sin  6 


1  -  e  cos  6 


£L^       •    -1 
+    sin 


o 


cos  6 


e  cos  6 


/gm! 


(4-22) 


and 


sNl-e     sin  6 


A 


1  -  e  cos  6 


+    sin 


A 


e  -  cos  8 


A 


1  -  e  cos  6 


A 


*f-*e 


GM, 


T  (4-23) 


Substituting  t„  =  t,  -  t     after  subtracting  (4-23)  from(4- 22)  we 
have 
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)\fl- 


sine 


o 


1-e  cos  G 


o 


sin  0. 
-  e  cos 


A 


+ 


sin 


-1 


/ 


cos  6 


e-  cos  6 


1-e  cos  0 


•sin 


A 


=  t,    7(-)b    n/GM 


U-e  cos  0 


(4-24) 


A 


where  GM 


0.  01716  if  a  is  given  in  AU.   t«  is  in  days  (sidereal) 


Values  for  a  large  number  of  trajectory  time  of  flights  were 
calculated  and  the  results  plotted  in  Fig.   4.10,     Fig.   4.2  illus- 
trates pertinent  factors  in  the  foregoing  equations  . 

4.12  Angular  Separation  of  Earth  and  Mars 

Having  found  the  time  of  flight  for  a  given  trajectory  we  may 
compute  the  angular  separation  of  Earth  and  Mars  at  the  instant 
of  firing  (6     -  0   ).     Since  Mars  is  in  motion  about  the  Sun(0    -  6   ) 

s         O  !3  OS 

is  a  quantity  which  must  be  known  in  order  to  fire  at  the  correct 
time  to  produce  an  intercept.     It  is  evident  from  Definition  Summary  1 
that 


°s    =  6A  +  ?.^   *f 


(4-25) 


6. 7 


It  may  be  shown    "  that 


\ir 


(1-e  cos  6  :. 


r'.2 


cr*     T      a_e2}s/2 


(4-26) 


=  constant  (1-e  cos  c^-*) 


(4-27) 


Since  6U^is  a  function  of  6U  pfwe  are  inexorably  engaged  with 
a  quantity  which  varies  not  only  with  real  time  but  also  with  ang- 
ular  position.     A  convenient  method  of  expressing  fy^v-^f  is  "co  plot 
Gr^tf  as  a  function  of  time;  then,    knowing  t-  and  choosing  t     we  can 
obtain  the  angle 


e    *  x  =  e    -  e  A 
©of        s       A 


(4-28) 
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and 


(So  -  V  ■    (6o  -  V  "  <9s  -  V  (4-29) 

The  aforementioned  plot  appears  as    Fig.     4,  11    .  •   To  assist 
in  the  logical  selection  of  t         an  overlay  (  Fig.    4.  12;  in  en- 
velope attached  to  back  cover)  has  been  constructed.     The  basis 
of  construction  and  description  of  use  is  given  in  Information 
Summary    4.1 
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Fig.   4.11     Mars  angular  position  measured  from  Earth  aphelion 
as  a  function  of  time 
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Fig.   4.12    (Overlay  #1) 


Earth  angular  position  measured  from  earth 
aphelion  as  a  function  of  time 


(In  pocket  on  back  cover) 
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INFORMATION  SUMMARY  4,  1 


Construction  and  Use  of  Overlay  #1. 

Assuming  the  "state  of  the  art"  of  interplanetary  travel  is 

such  that  an  Earth-Mars  trip  is  planned  to  occur  about  1980, 

2    54  o 

we  note   J  ""     that  on  January  1,   197  9  Earth  is  177     ahead  of  her 

aphelion  (see  insert  on    Fig       4.12    ).     Mars  is  2     behind 

Earth's  aphelion,     In  other  words  Mars  is  181     ahead  of  the 

Earth o     During  the  succeeding  1,135  sidereal  years  Earth 

overtakes  Mars.    and*they  are  in  conjunction  227     past  Earth 

aphelion  (about  February  17     1980).     Hence  from  January  1,  1979 

to  February  17.    1980  (6     ~  0   )  takes    on  all  values  from  0  to  7r, 
J  o         s 

With  all  pertinent  parameters  of  both  Earth  and  Mars  orbits 

54 
known  ;,      Fig.   4.  II  and  4.12  were  drawn  so  that  information  for 

the  period  in  question  would  be  easily  accessible.     The  as- 
sumption has  been  implicitly  made  herein  that  Mars  aphelion 

precesses  only  slightly  about  the  Sun  (a  few  seconds  of  arc  in 

54 
about  10, 000  years). 


(a)  Basis  for  Construction  of  Overlay 

The  angular  position  of  Earth  is  indicated  for  January  1,1955 
and  will  repeat  every  sidereal  year       Mars  angular  position  is 
marked  on  graph  4.11    for  January  1,    1955;    but  its  period  is  not 
an  even  multiple  of  years,   hence  its  angular  position  for  any 
given  time  must  be  found  either  from  an  Ephemeris  or  by  the 
following  method: 


(Page  1  of  3] 
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Divide  the  difference  in  time  between  January  1,    1955  and 

the  time  for  which  Mars  angular  position  is  desired  by  the 

number  1.  88  (Mars  period  in  sidereal  years).     Slide  the  mark 

for  January  1,    1955  on  the  overlay  to  the  right  of  the  mark  for 

January  1,    1955  on  the  underlying  graph    an  amount  equal  in 

sidereal  years  to  the  remainder  in  the  foregoing  division.     The 

real  time  scales  have  now  been  superimposed  so  that  the 

angular  positions  of  Earth  and  Mars  are  coincident  on  the  time 

scale.     The  values  of   0     and  0     may  be  picked  off  the  curves 

o  s         J         * 

and  ordinate  scale.     If  a  time  other  than  January  1,   for  a  given 
year  is  desired,    make  auxilliary  pencil  marks  for  Earth  position 
and  use  that  rather  than  January  1,    as  a  reference. 


(b)  To  find  (0     -  0   )  for  a  given  t 
o         s _ o 

During  the  period  of  interest  all  values  of  (0     -  0   )  from 
&  ^  OS 

0  to   7r  occur,   hence  set  the  overlay  as  before  and  read  t     corre- 
sponding to  the  given   0     and  0     as  a  decimal  fraction  of  a  year 

o  s 

after  the  chosen  reference  time. 


(c)  To  find  t     for  a  given  0     and  0 
o ° o s_ 

This  problem  is  somewhat  more  complex.     We  are  primar- 
ily    concerned  with  a  real  time  firing  time  in  late  197  9.     This 
is  true  because  it  is  about  this  time  that  (0     -  0   )  takes  on 
values  in  which  v/e  are  interested  (note  Section  4  12  and  Fig. 
4.11        and    4.12     )-    Fortuitously;  Mars'  angular  rate  in  late 
1979  is  approximately  constant  at  0.47  degrees  per  sidereal 
day.     If  we  multiply  0.  47  times  the  tf  this  determines  a  rough 
value  of  (0     -0A).     Subtracting  (0     -  0A)  from  (0     -  0A)  pro- 

S  J\  S  J\  O  x\ 

vides  (0     -■  0   )  for  which  a  t     can  be  found  as  explained  in  part 
os  o  r  tr 

(c).     By  nothing  the  angular  positions  of  Earth  and  Mars  and  by 

moving  a  straight-edge  up  or  down  the  superimposed  curves 

one  may  obtain  a  more  accurate  value  of  (0    -  0   )  for  the  given 

data. 

(Page  2  of  3)  
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(d)  To  find  t     for  a  given  t-,    6     and  6  . 
o  &  V      o  A 


January  1,   1955  was  chosen  as  reference  time  because 
accurate  values  were  available  to  the   authors . 

The  method  described  herein  is  a  rough  graphical  method 
of  finding  required  quantities  for  preliminary  work.     Actual 
calculation  would  of  necessity  be  made  by  computer. 


(e)  Notes 

(Page  3  of  3) 


4.13    VeA  and  Vq?  .  as  functions  of  (0     -  0g) 


Since  we  now  may  find  (6     -  6   )  for  any  given  trajectory, 
it  appears  pertinent  to  note  any  relations  which  may  exist 
between  this  angle  and  the  velocity  required  to  fly  the  trajecto- 
ry.    From  an  inspection  of  the  basic  equations  in  Section  4.8 
it  is  obvious  that  the  relationship  is  quite  complex.     For  this 
reason  V^A  and  V^.  were  plotted  graphically  for  a  large 

number  of  trajectories  in  a  manner  analogous  to  that  of 

51 
Preston-Thomas     .     The  minimum  V^  A  was  found  to  be  the 

envelope  of  the  curves  as  6     took  on  various  values.  (See  Fig. 

4. 13     )     The  minimum  V^A  was  found  to  be  identically  the 

curve  for  6 «  =  0    .    (See    Fig.    4.14   )  This  is  not  inconsistent 

since  the  first  case  represented  minimum  velocities  for  a 

non-maximum  radius  vector  of  an  ellipse;  whereas  the  second 

case  represented  the  minimum  velocities  that  exist  at  aphelion 

(maximum  radius)  of  an  ellipse.     In  both  cases  the  Hohmann 

Ellipse  (6     =  180   ;  G  A-  0   )  provided  the  absolute  minimum 
o  J\ 

velocities. 


4.14    The  Rocket  Plane  (A  plane)  Concept 

Let  us  define  the  plane  containing  Rq/»  and  V^,     at  t  =  t 
as  the  rocket  plane  (A  plane).     If  the  rocket  were  not  perturbed 
during  its  flight  it  could  remain  in  this  A  plane  (this  being  a 
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Fig.   4.13     Rocket  velocity  magnitude  with  respect  to  Earth  as  a  function 

of  6     and  6-6 
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44 
property  of  the  motion  of  a  body  in  a  central  force  field     ). 

The  A  plane  need  not  necessarily,    ana  generally  will  not,  lie  m 
the  Earth  ecliptic  plane  since  Mars  passes  through  the  Earth 
ecliptic  plane  only  twice  every  1.  88  sidereal  years.    Unfortu- 
nately,  for  calculation  purposes,   the  orientation  of  Vq  a  in 
Solar  System  space  depends  on  the  real  time  firing  time  and 
time  of  flight.     Hence  every  possible  trajectory  requires  a 
complete  individual  solution.     The  labor  involved  is  so  moun- 
tainous that  solution  by  high-speed  digital  computer  is  manda- 
tory; and  even  these  solutions  must  be  based  on  a  selected 
future  firing  time 

If  we  can  show  that  the  displacement  of  the  rocket  normal 
to  the  A  plane  is  negligibly  small  we  may  be  justified  in  saying 
that  for  practical  pusposes  the  rocket  remains  essentially  in 
the  rocket  plane  and  our  problem  reduces  to  two  dimensions  , 
Neglecting  the  possible  perturbations  of  bodies  such  as  meteors, 
distant  planets,    asteroids,    etc.,   the  bodies  causing  displacement 
of  the  rocket  from  the  A  plane  are  Earth,    Moon,    and  Mars  .     In 
an  effort  to  show  the  minor  nature  of  the  aforementioned  dis- 
placements calculations  were  made  on  a  single  "typical" 

trajectory.     This  trajectory  was  chosen  because  V.        =10,  000 
j  j  j  j  inc 

ft /sec,    which  is  a  velocity  probably  within  the  capabilities 
achievable  in  the  near  future      Since  we  are  concerned  with 
the  period  from  late  1979  to  early  1930  the  firing  time  is  es- 
tablished by  (0     -  d   )  to  be  about  December  6,    1979      Note 
J       o         s 

Figs.  4,15  and  4.18. 

Referring  to  Derivation  Summary  4   1  the  cumulative  dis- 
placement of  the  rocket  from  the  A  plane  was  found  by  ac- 
cumulating the  values  determined  from  Equation    11  over  the 
trajectory  traveled  by  an  unperturbed  rocket(see     Fig    4.  19  ). 
Certain  errors  are  inherent  in  the  calculations  but  will  be  seen 
to  be  negligibly  small.     The  arc  length  of  the  trajectory 
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ROCKET  ELLIPSE 
APHELION 


DATE  OF  FIRING 
APPROX.  DEC.  6,1979 


MARS 

ASCENDING 

NODE 


MARS 
PERIHELION 


Fig.   4.15      Angular  quantities  of  a  representative 
trajectory  discussed  in  section  4.14 
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Fig.   4.  16      A  typical  trajectory  showing  the  effect  of  tilt  to  account  for 
non-co-planarity  of  Earth  and  Mars  motion 
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AN   IMAGINARY  CYLINDER 

RADIUS   EQUAL  TO  rrf  ,  Z 

PERPENDICULAR  TO   E 


*  e^pse 


%*$■ 


set 


*■      -'Aplane 
measured  from  rocket 


plane 


Fig.   4.17      Assisting  supplementary  view  to  figure  4.  16 
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DERIVATION  OF  THE  EQUATION  FOR  THE  DETERMINATION 
OF  THE  DISPLACEMENT  OF  THE  ROCKET  NORMAL  TO  THE 

"ROCKET  PLANE"         "» 


Four  relations  are  available  from  consideration  of  the  geo- 
metrical properties  of  the  trajectories  or  orbits.     Two  assumptions 
are  made  which  greatly  simplify  the  solution.     The  error  introduced 
by  these  assumptions  is  on  the  "high"  side  and  will  be  shown  to  be 
negligible  (Fig.    4.16  refers). 

sin  <5 


(i         )A 

max  A 


icf=  'We?31" ik  +  eo  -  6d  > 

z     =  R~  A  sin  iA 
A         CA  A 


zd^  rd sin  Hi 


(i) 

(2) 
(3) 
(4) 


From  (3)  and  (4) 


i.    =  sin 
A 


si 


R 


sm  l 


CA 


d 


(5) 


(a)    Basic  Equations 


To  aid  in  clarification  of  the  development  that  follows  and 
the  assumptions  made,   the  following  list  of  symbols  is  presented: 

(g      ,.  ).  Magnitude  of  that  component  of  G.  normal  to 

J  J 

the  "rocket  plane". 


DERIVATION  SUMMARY  4.1 
(Page  1  of  4) 
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i.  Inclination  of  indicated  quantity  with  reference 

plane,   or  with  ecliptic  plane  if  no  reference 
is  given-positive  direction  in  the  conventional 
sense. 


(i         ).      Maximum  inclination  as  defined  above 

max  j 


z.  Displacement  of  indicated  quantity  normal  to 

indicated  reference  or  to  ecliptic  plane  if  no 
reference  is  given. 


A  Angle  as  shown  in  Fig.  4.16 

£"  Eccentric  anomaly  as  shown  in  Fig.  ±.  16 

6  Increment  (Magnitude  as  defined). 

(b)  Definitions 

Referring  to  Fig.  4.  16  and  4. 17 the  two  assumptions  made  are- 
(  1  )  that  the  "rocket  ellipse"  is  almost  circular  about  its 
center  (C). 

(2  )  that  the  line  CC  in  Fig.  4.  16     is  negligibly  small. 

(c)  Assumptions 

To  derive  the  required  equation  we  solve  (      3         ),   then 
(      4      )  for  the  final  condition  &■<=  6..     We  then  find  i.  from 

(     5     )  and  (i         )A  from  (1  ).     Having  discovered  (i         )A 

max  A  &  max  A 

for  a  given  trajectory  we  are  now  able  to  write  i.  as  a  function 
of  Gtf  (Equation    1  ) . 

Any  [  z     ]  A  i  may  be  found  thusly: 

DERIVATION  SUMMARY  4.1 
(  Page  2  of  4  ) 
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From  Fig,  4„  17       ; 


fZ.    ,        ]        *  -    \  R~  A  s  in  i  A 1        j, 
1     AplaneJ  at^J      L     CA  AJ  at^J 


and 


obviously 


zo*  =  ir<5sin  tfUt 


6 


(6) 


(7) 


[  Zdl  Aplane=  [  ZAplane]  at^  "  Z<3 


(8) 


now 


fid) 


Aplane 

measured  from 
rocket 


sin 


-1   iZd  J  Aplane 


T&    r, 


(9) 


then 


and 


[Sorthlc?    =      ^at(ld-   rA> 


[Z^l 


Aplane 


\r&    rA 


r    6  ZA  1    A    .  (due  to  Mars) 

L  AJ  Aplane 


[  gorth]C?  (6t)2 


(10) 


[gorth^(6t)2x7"15xl°5 


(ID 


if  6  t  is  in  sidereal  days 

g     ,.     is  in  ft/sec 
&orth 

6  Z  A  is  in  miles 

DERIVATION  SUMMARY  4.  1 
(Page  3  of  4) 
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A  similar  derivation  may  be  made  for  the  cis -lunar 
Barycenter. 

(d)    Derivation 


DERIVATION  SUMMARY  4.1 
(Page  4  of  4) 


traveled  by  the  rocket  is  approximately  354,  000,  000  miles;  the 
cumulative  displacement  from  the  A  plane  is  approximately  5,  500 
miles:  a  ratio  of  about  1.  55  to  100,  000,    certainly  negligible. 
Further,   the  calculations  show  that  the  displacement  due  to  the 
mass  attraction  of  Mars  is  only  33  miles,   the  remainder  being 
the  effect  of  the  Earth-Moon  mass  attraction.     About  3,250 
miles  displacement  occurs  in  the  first  five  days  of  travel;  hence, 
if  a  computer  program  were  established  to  counteract  this,   the 
displacement  could  be  reduced  even  further. 

The  obvious  advantages  of  reducing  the  problem  to  two 
dimensions  are  many -fold.     These  include  (a)  reduction  in 
weight  and  complexity  of  the  guidance  mechanism;  (b)  simpli- 
fication of  component  construction;  (c)  savings  in  fuel;  (d)  im- 
provement of  payload  ratio;  (e)  saving  in  material,    labor,    and 
cost;  (f)  contraction  of  production  time-tables;  and  many  others. 


78 


CHAPTER  5 

DEPARTURE  AND  ARRIVAL 
ORBITS 


5,1      Introduction 

The  capability  of  a  rocket  to  achieve  required  velocities  de- 
pends so  critically  upon  the  energy  carried  within  itself  that 
every  effort  must  be  made  to  minimize  these  energy  requirements 
This  section  investigates  one  aspect  of  this  minimization  --  the 
selection  of  departure  and  arrival  orbits  about  the  respective 
planets  to  minimize  kinetic  energy,    hence  velocity,    that  the 
rocket  must  provide  itself  in  order  to  fly  a  given  trajectory. 
Fig.   5.1  illustrates  the  relationships  of  pertinent  velocities . 

For  complete  definitions  see  Definition  Summary  1 


inc      ®s 


Fig,   5.1      Pertinent  Velocity  Relations 
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5.  2     Optimum  Orbit  Radii 

In  order  to  fly  any  given  trajectory,    or  ellipse,   a  velocity  of 
the  rocket  with  respect  to  the  Sun  (Vq  .)  must  somehow  be  attained. 
Referring  to  Section  4.9,   we  note  that  this  is  the  velocity  vector 
required  for  an  intercept  of  the  orbit  of  Mars,   treating  both  Earth 
and  Mars  initially  as  massless  points.     However,   this  is  obviously 
not  the  case.     Referring  to     Fig.    5.  2  we  can  say,   figuratively, 
that  we  must  climb  up  the  side  of  the  Earth  potential  "well",   over 
the  intervening  "hump",   then  down  into  Mars  potential  "well. 
Jupiter's  potential  "well"  is  shown  for  comparison.     From  Fig.  5.  1 
we  note  that  the  rocket  already  possesses  a  certain  velocity  with 
respect  to  the  Sun  by  virtue  of  the  Earth's  orbital  motion,   thus  the 
rocket  must  be  provided  the  velocity,  Vm  ..     In  order  to  "over- 
come" the  effect  of  the  earth  gravity  field,   the  rocket  must  be 
provided  a  greater  velocity,   V,^Av    .     The  rocket  also  possesses 
a  velocity,   V^     ,   by  virtue  of  being  in  orbit  about  the  earth, 
therefore  it  remains  for  the  rocket  to  provide  by  its  own  efforts 
an  increment  of  velocity  V.      .     Initally  Vq  a  will  be  along  the 
dotted  line  labeled  in  Fig.    5.1,   however,   the  earth  gravity  effect 
decreases  rapidly  so  that  VX   .    lies  almost  identically  on  the  solid 
line.     In  addition  the  ratio  \rn\    /R^va    is  of  the  order  10      ,   hence 
the  displacement  of  the  rocket  from  the  earth  rgr\     is  negligible 
for  purposes  of  calculating  initial  velocity  requirements. 

Regardless  of  validity,    if  Vq  a  is  known  we  may  obtain 

relations  between  V.       and  rjs    .     This  has  been  done  in  Derivation 

inc  Ws 

Summary  5.1.  Although  the  relations  have  been  developed  for  the 
planet  Earth  the  derivation  is  perfectly  general  and  applies  to  any 
body;  with  different  constants. 

Parts  (b)  and  (c)  of  Derivation  Summary  5.1  develop  and 
prove  the  optimal  relationship  of  V.        and  rg\    .    Part  (d)  expresses 
the  relation  non-dimensionally. 

Pertinent  ranges  of  V.        vs  rm     have  been  plotted  in  Fig. 
b  inc  w  s  to 

5.  3  and  the  non-dimensional  function  is  shown  over  an  applicable 
range  in    Fig.     5.4. 
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IN    FT/SEC x 10 


Fig.  5.3      Optimal  relationships  between  V.       and  r.     for 

J 

Earth  and  Mars 
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is  the  frequency  of  occurence  of  the  rather  special  Earth-Mars 
configurations  required  bythe  Hohmann  ellipse.     This,   among 
other  considerations  has  been  taken  up  in  preceding  sections. 


H4 


5.  3      Use  of  Optimum  Orbit  Concept 

Almost  paradoxially,   the  optimum  orbit  concept  is  of 
particular  use  in  the  selection  of  non=optimum  orbits  or  tra- 
jectories.    As  an  example,    assume  that  a  rocket  has  been 
developed  with  a  capacity  to  provide  itself  with  a  V.       of 
approximately  10,  000  ft /sec       The  optimal  orbit  radius  corre- 
sponding to  10,  000  ft/sec,    as  an  optimal  velocity,    from  graph  5    3 
is  approximately  27,  500  miles  above  Earth  center.    If  27,  500 
miles  is  unacceptable  due  to  other  considerations  we  can  reduce 
(from  graph    5.4    )  this  radius  by  42.  5%  simply  by  allowing  a  2% 

increase  in  the  required  velocity.     In  other  words,   for  V. 

n  J  inc 

=  10,  200  ft /sec    r  ~     =15,  800  miles,   or  approximately  12,  000 

miles  above  the  Earth's  surface.     It  must  be  borne  in  mind, 

however,   that  these  relations  exist  only  for  that  trajectory  for 

which  (V.      )      .  =  10,000  ft/sec 
inc  opt 

Applying  the  foregoing  reasoning  to  the  "Hohmann"  type 

ellipse  -(V.      )         =6,  818  ft/sec;  (r  m   )      .  =  57,  300  miles  -  we 
^  inc  opt         '  (Bs  opt 

may  reduce  r~     to  33   000  miles  by  increasing  V.       to  6,  936 

J  0s  J  &      inc 

ft/sec.     Hence     if  we  can  provide  the  energy  to  increase  (V.      ) 

r  &J  inc  opt 

by  118  ft /sec  we  not  only  gain  the  obvious  advantages  of  a  lesser 
r  q     but  also  the  flexibility  inherent  in  a  greater  freedom  of 
choice  in  the  angular  positions  of  Earth  and  Mars.     This  has  been 
discussed  more  fully  in  Chapter    4 


5   4      The  Least  Orbit  Concept 

If  the  problem  is  posed  that  V.        =  7    500  ft /sec.  ,  what  is 

^  r  inc  ' 

the  least  orbit  the  rocket  may  depart  from  and  still  travel  a 

trajectory  to  arrive  at  Mars  ?     This  obviously  is  not  an  optimum 

orbit  problem  as  previously  posed      For  (V.      )         =  7,  500  ft/sec 

(r  -,   )         =  47,000  miles  (graph  5.3).     For  (V.      )         =6,818 
®s  opt  &     ^  inc  opt 

ft /sec  (the  minimum  "Hohmann"  ellipse  velocity)  the  ratio  V. 

^  J  mc 

/ (V .      )      ,  =  1.1  and  from  graph  5.  4     the  ratio  r  ,~     / (r  ^    ) 

inc  opt  &      ^  0S'       ©s  opt 

=    0.  327  and  r  ~      =  18,  750  miles       The  limiting  restriction  herein 
^p  s 
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DEPARTURE  AND  ARRIVAL 
ORBITS 


Optimum  Planetary  Orbits 

In  order  to  "escape"  any  bodies  mass  attraction  a  vehicle 
must  attain  kenetic  energy  at  least  equal  to  the  potential  energy 
inherent  in  the  field  or  fields  to  be  "overcome";  hence 

KE  =  mA  (  |G.A    dR  (1A) 

A,)  I     jA 

or 

r.oo      GM 
^      I        3   dR  (IB) 


mA        '  R      R2 
Following  the  convention  Q  =     |q|    we  integrate 

and  ^T^tH  <2> 

111  A  K  •   A 

A         jA 


but  —  =  1/2  V2.  (3) 

mA  Je 


XA 


/2GM. 
hence  V.     ^-,/tt- ^  (4) 

Je     vHiA 

-v 

c  o 

It  has  been  shown       that  a  rocket  in  circular  orbit  about  a 
body  possesses  a  so-called  circular  velocity. 


DERIVATION  SUMMARY  5.1 
(Page  1  of  5) 
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/GM. 
V      -  V^-J  (5) 

hence 


V.     =  nHTv.  (6) 

]e  js 

From  Fig.       in  the  case  of  the  earth  and  a  rocket 

V/mA*     =  V.        +  V~  (7) 

(®A)g         inc        ©s 

Further,   by  equating  the  kinetic  energy  requirements  for 

the  rocket  we  have  . 

1/2  mAV     Al/2mAV       2  =  l/2mAV2,^M 
A0    A  A©  e  A®  (®A)g 

(8) 

where  the  first  term  represents  the  kinetic  energy  requirement 
to  achieve  the  velocity  "vV^A;  the  second  term  preresents  the 
kinetic  energy  requirement  to  "escape"  the  earth  gravity  field; 
and  the  third  term  describes  the  sum  of  the  first  and  second. 

From  (8)  we  have 


V(©A)g=  ^AtV*e  (9> 

Hence,    substituting  in  (7)  and  rearranging 


V         =  Vvl  "  +  2v£  -  Vm  (10) 


inc  "    vv0  A  x     v0  s  v0s 

and  if  initially  R     .   =  r     „  then 


©  A    r©s       res 


(a)  Development  of  V        in  terms  of  rocket  orbital  radius 

r  inc 


DERIVATION  SUMMARY  5.1 
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In  order  that  r^     be  the  optimum  radius  on  a  kinetic  energy 

basis  V.       must  be  a  minimum  for  r/T>     =(rcu   )      L  which  is  the 
mc  Ws        a?s  opt 


equivalent  of  the  conditions  <•> 

6V. 

(1)  7-^    =  0  (12) 

6r0s 

62V. 

(2)  2mc-  >    0  (IS) 

6r0  s 

To  satisfy  (1)  we  take  the  first  partial  derivative  of  (11) 
and  equate  to  zero. 

GM0  GM0 


6V. 

mc 

2 

2 
r0  s 

6r0s 

1       ™e"'- 

V V^  a+ 

0A    r0s 

0S 

0  (14) 


Solving  for  r0g  =  (ifc^ 
2 


2GM  2    '       _         /^% 

^Wopt+V®A  <Wopt 


2GM$ 
's'opt 


and  as  corollaries 


(15) 


tt?s  oot  „2 

V0  A 


VmA    =-n/T«.V/k    )     ,  (17A)         (V.      )         =   / —        (17B) 

tDA  n?s  opt  mc  opt     J  rm 

(b)  Derivation  of  (r^   )  optimum 

DERIVATION  SUMMARY  5.1 
(Page  3  of  5) 


To  satisfy  condition  (2)  a.bove  we  take  the  second  partial 

derivative  of  V.       with  respect  to  r^     and  evaluate  at  rm^, 
inc  ^  ws  <£>s 

=  (rg   )        to  show  it  is  indeed  positive , 


62  V. 


inc 


x     2 

6  re 


G% 


©  s 


GM 


-    -2 


© 


GM 


© 


2 


©     s 


2GM0 


res 


+  V, 


3/2        /GM0\3/2 


©  A 


© 


I 


(18) 


+ 


1 


+    V/T\ 


A 


ngnvge 

res 


Combining 


GM0V 


62V 


inc 


re 


6r; 


2GM6>     v2       U/2 

^7+ VH 


GM^ 


"372" 


re; 


(19) 


!GM0 


+ 


re 


2GM©    tr2 

+  V©   A 


r©s 


Evaluating  (19)  for  rm 


2GM0 


v4a 


GM® 

3 

0  s 

CM® 
6>s 


or  V> 


2GM0 


©    A 
DERIVATION  SUMMARY  5,1 
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62V. 

inc 

6r©  s 

1      v4                        I      V4 
4     V©   A                   8      V©  A 

2^   resV^A     2^re  s    V0A 

v 

(20) 

V2                                   ^-V2 
V©   A                           2  V©   A 

^   r©s    V©A     n^4  sV©A 

V©A    .      1                 1         . 
2          K    2n/2"           8n/2"    ' 
r©  s 

(21) 

3V©A 

(22) 

8^2   r|g 

which  is  positive  for  positive  values  of  the  radical 

(c)  Proof  of 

Minimality  of  (r^)^ 

Dividing 

(11)  by  (17B)  and  noting  (17A)    we  have 

I 

V.                     /     2                            /  1       ' 
mc            _/                     0           1 

(23) 

^nc^pt      /r©s                         /   r©s 

<s/(r©fe}opt             «/(lWqpt 

(d)  Non-dimensional  Relations 

DERIVATION  SUMMARY  5.1 
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CHAPTER  6 

THRUST  AND  ENERGY  CONSIDERATIONS 

A,    THRUST  REQUIREMENTS 

6.  1    Introduction 

It  is  of  interest  to  examine  the  thrust  requirements  and 
pay  load  ratios  for  a  rocket  departing  circular  orbit  about  earth. 
Such   figures  appear  to  be  more  or  less  well  known  for  the 
journey  from  earth's  surface  into  orbit  (as  for  earth  satellites). 
However,    the  orbit-to-orbit  requirements  are  not  commonly 
available,   although  easily  calculable.     This  information  is  con- 
sistent with  the  optimal  orbit  concept  developed  previously  and 
verifies    the  validity  of  such  an  approach  for  the  conditions 
assumed 

Although  the  data  have  been  based  upon  the  departure  and 
arrival  orbit  concept,   an  extension  to  surface  departure  can 
easily  be  made  by  extrapolation  to  lower  altitudes.     The  penalty 
of  departing  a  lower  altitude  will  be  evident  after  an  investigation 
of  this  section 

6.  2  The  Optimal  Payload  Ratio 


The  change  in  velocity  necessary  to  affect  transfer  from 
orbital  motion  about  earth  to  departure  path  has  been  designated, 
V.      .     This  change  in  velocity  is  accomplished  by  applying  thrust 
for  a  controlled  period  of  time,     If  the  power-on  period  is  suf- 
ficiently brief    the  maneuver  trajectory  can  be  neglected  and  the 
transition  from  circular  orbit  to  departure  path  thought  to  occur 
instantaneously  at  tangency.     The  justification  of  this  assumption 
will  be  discussed  shortly  after  the  development  of  the  equations. 

The  commonly   used  expressions  for  an  idealized  single 
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step  thrust  study  of  this  nature  are  derived  from  the  equations  of 
motion  written  along  the  orbit  flight  path. 

m~dt~  =  T  -  D  -  Mg  sin  6 

In  circular  orbit,   however,   for  the  above  assumptions,   D 
is  presumed  negligible  and  0=0. 

m^7=  T  =    M     -   -EL       -S   for   0<t<t 
dt  o        t,  dt  -    -  b 

b 

where    M     =  initial  gross  mass 
o  6 

M     =  mass  of  propellant 
t,      =  time  of  burning 

A  constant  thrust  device  is  assumed  of  specific  impulse,   I. 


M 

T  - 

dv 

dt  " 

M 

tpig 

xb 

M    -M 

b 

t — =—     where 

l-o  t 

Mp 

M 
o 

Integrating  between  t  =  0  and  t  =  t. 


V     -  V     =  V.       =   -  Ig  In  (1  -  6)=   -  Ig  In  (PR) 

t,  o  inc  s  »  s 

b 


PR    =  Payload  Ratio  = 


M     -  M 
2         p 
M 
o 


It  is  to  be  remembered  that  this,   V.  ,    applies  only  to 

inc@  '      *^  J 

the  departure  maneuver,   as  does  PR(g  ■     For  entry  into  circular 
orbit  about  Mars,   a  similar  expression  is  obtained.     The  total 
payload  ratio (  PR.)  for  an  orbit-to-orbit  maneuver,    ignoring 
corrective  thrust  requirements,    is 


PRt  =  PR#  PR  a 

These  equations  have  been  expressed  in  tabular  and 
graphical  form  for  various  values  of  specific  impulse,   which 
has  been  chosen  as  the  propellant  parameter.    The  subscripts 
in  the  headings  of  Table  6. 1     refer  to  the  three  assumed  values 
of  specific  impulse.     Table  6,  If  or  example,   lists  Payload  ratios 
for  various  sample  trajectories  developed  by  the  methods  of 
chapter    4     .     The  importance  of  improved  propellant  performance 
is  noted.     Although  the  corrective  thrust  requirements  have  not 
been  included,    reference  47 


e    -  e 

o         s 

*f 

Vinc© 

PRe250 

PR©300 

PR®350 

55.8° 

201 

7391 

.022 

.042 

.063 

58.2°* 

259.6 

6818 

.065 

.101 

.139 

21.7° 

145.8 

23,  696 

.004 

.011 

.021 

16.2 

125.5 

24,  529 

.002 

.005 

.011 

20.1 

181  .0 

19, 045 

.009 

.021 

.037 

27.0 

134.5 

16, 763 

.003 

.008 

.017 

7.0 

241.5 

15,181 

.018 

.036 

.057 

33.2 

153.1 

13,055 

.008 

.017 

.030 

39.3 

204.1 

11,  932 

.031 

.054 

.082 

43.0 

228.2 

9358 

.045 

.075 

.108 

48.1 

185.4 

8151 

.019 

.037 

.060 

45.3 

255.4 

7290 

.060 

.095 

.134 

38.7* 

171.5 

10,146 

.005 

,011 

.021 

TABLE  6.1 
SOME  REPRESENTATIVE  TRAJECTORIES  WITH  CORRESPONDING 

PAYLOAD  REQUIREMENTS 
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has  indicated  that  a  suitable  payload  ratio  for  corrective  thrust 
requirements  is  .79.     If  the  above  figures  are  altered  by  this 
multiplier,  it  is  seen  that  the  most  optimistic  Payload  Ratio  to 
be  hoped  for  is  ,11.     The  calculations  of  Table   6,1  presume  de- 
parture and  arrival  orbits  at  "optimal"  altitude  as  defined 
earlier.     Since  these  altitudes  inevitably  vary  for  various 
trajectories,   little  correlation  is  to  be  expected  between  Payload 
Ratio    (PRrp)  and  angular  separation  at  instant  of  fire  (6     -  6   ) 
other  than  to  indicate  a  lower  bound  for  certain  specified    para- 
meters , 

6.  3  An  Interpretation  of  "Optimum"  Altitude 

At  first  glance,    it  appears  that  the  higher  the  optimal 
departure  and  arrival  orbit  altitudes,   the  more  favorable  the 
total  Payload  Ratio.     To  invalidate  this  idea  as  well  as  to  demon- 
strate the  advisability  of  employing  optimal  altitude  two    tra- 
jectories have  been  selected  for  study:  one  is  the  familiar  Hohmann 
semi-ellipse  trajectory:  the  other  merely  a  representative  one 
selected  for  contrast.     Both  are  defined  by  the  *  in  Table  I. 

If  the  PR(g     and  PR^f    for  these  two  trajectories  are 

calculated  as  functions  of  altitude  (including  optimal  altitude)  on 

a  non-dimensional  basis,   as  in  Fig. 6.  land  6.2    it  can  be  seen 

that  a  peak  payload  ratio  is  achieved  at  precisely (r  )/~)prr.=  A  slight 

fall -off  occurs  at  values  offr/  >  1-      Thus,   the  choice  of 

s  OPT 

optimum  altitudes,    in  the  sense  defined,    insures  a  maximum 

total  payload  ratio  for  any  given  specification  of  V.      ^  and 

mc  tt? 

V.      ^r .     The  great  saving  in  Payload  Ratio  is  noted  when  one 

departs  from  circular  orbit  altitude  for   r^-\  =  .5. 

°VOPT 

6.  4  The  Assumption  of  Impulsive  Thrust 

A  justification  of  the  assumption  of  impulsive  thrust  is 
dependent  upon  a  comparison  of  burning  time  and  period  of  orbit. 
Obviously,   for  a  low  thrust  device  of  considerable  duration,   the 


PR©  .3 
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Fig.  6.1     Payload  ratios  VS.  orbit  altitude  for  the  Hohmann  ellipse 
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PR*  .3 
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^0-6^=  120° 
1 

1 
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Fig.  6.2    Payload  ratios  vs.  orbit  altitude  for  a  representative  trajectory 
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Fig.  6.3     Orbit  period  as  a  function  of  altitude  above  planet  center 


97 


assumption  is  not  valid. 

The  period  of  orbit  for  both  Mars  and  Earth  is  plotted  in 
Fig,  6.  3  as  a  function  of  distance  from  Earth  and  Mars  center. 
If  we  assume  that  instrumentation  is  capable  of  withstanding  an 
axial  load  of  lOg  upon  departure  from  Earth  orbit  the  minimum 
acceptable  burning  time  can  be  determined.         Careful  design 
can  then  be  employed  to  produce  a  favorable  burning  time.     The 
minimum  burning  times  as  functions  of  orbital  period  for  optimum 
Earth  and  Mars  orbit  are  listed  in  Table  6.  2      for  the  trajectories 
previously  considered:    A  specific  impulse  of  350  seconds  is 
assumed. 


r 
s© 

T    (min) 

sec. 

48,  772 

3550 

33 

57, 295* 

4650 

29 

4,  744 

115 

252 

4,  424 

275 

7,  344 

220 

154 

9,  481 

320 

120 

11,  558 

400 

100 

15,  630 

680 

76 

18,  709 

880 

66 

30,  416 

1900 

45 

40, 099 

2750 

37 

51, 782 

4000 

32 

25, 877* 

1420 

51 

Ttf 

*b 

3353 

210 

2  54 

7792 

735 

102 

4571 

330 

171 

2585 

142 

334 

5471 

430 

142 

2127 

108 

438 

6315 

530 

123 

2607 

145 

334 

6918 

620 

113 

7404 

680 

108 

7689 

700 

247 

3261 

200 

102 

1631 

71 

640 

TABLE  6.  2 
A  COMPARISON  OF  ORBITAL  PERIOD  AND  BURNING  TIME 
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B.  ENERGY  CONSIDERATIONS 


6 .  5    Introduction 

It  would  seem  appropriate  at  this  stage  to  mention  opti- 
mization of  the  trajectory.     The  term,   optimization,   however, 
invariably  leads  to  ambiguity  and  requires  explanation.     Do  we 
desire  optimization  with  regard  to  time,   fuel  expenditure,   or 
initial  velocity?--  or  possibly  cost?    Previously,   an  attempt  was 
made  to  minimize  the  initial  velocity  requirements,   but  it  is 
clear  that  this  does  not  necessarily  represent  the  optimum  as 
far  as  other  variables  are  concerned. 

An  overriding  consideration  during  early  interplanetary 
flights  will  undoubtedly  be  minimization  of  fuel  expenditure.     By 
minimization  of  fuel  expenditure,   we  imply  a  determination  of 
not  only  the  minimum  duration  of  powered  flight  but  also  the 
direction  in  which  this  power  or  thrust  is  applied.     Accordingly, 
these  are   all  considerations  in  the  minimization  problem.     It  is 
the  intent  of  this  section  to  describe  the  problem  in  general  terms 

and  summarize,    without  proof,    the  findings  of  other  authors  who 

L  1        -xl.  o-U.  ,  u-  34^  35,  36>  37 

have  written  at  length  on  this  particular  subject. 

6.6  General  Theory 


Let  us  first  consider  the  problem  of  optimum  transfer  (with 
respect  to  fuel  expenditure)  between  two  terminals  lying  in  a 
plane,     If  we  consider  an  inertial  orthogonal  coordinate  system, 
the  equations  of  motion  for  the  rocket  are 

dM  (6t1) 

dt 

*M  (6.2) 

dt 
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•  • 

X 

+ 

G      = 

X 

'cosaI 

•  • 

Y 

-!- 

G     - 

y 

cosp  M 

Z    +    G 


-cos  7 


_c_    dM 
M     dt 


(6-3) 


where  G   ,   G   ,   and  G     are  the  components  of  the  resultant  at- 
tractive  force  per  unit  mass  and  a^    p  and     y    are  the  angles 
between  rocket  thrust  direction  and  coordinate  axes.     The  ex- 
pression -jT^r   -r—   is  the  well  known  simplified  term  for  rocket 
thrust  and  is  defined  in  the  section  describing  thrust   consider- 
ations ,     c  is  the  effective  exhaust  velocity  and  M  is  the  mass  of 
the  rocket  at  any  time  t. 

G   ,   G   ,   G     depend  explicitly  on  x,   y,    z,   and  t  only, 
x        y         z 

Squaring  and  adding  6-1,    6-2,   and  6-3. 


[  X  +  Gx]  \[  Y+Gy]  2+  [  Z+GJ  2=  [  £  ffi]  2[  cos2  a+  cos2^  cos2  y] 

(6-4) 

1/2 


Integrating, 

c  log 

&e  M 


M 


o 


1 


[X+  G   1    +[Y  +G   ]    +[  Z  +  G  ] 


dt 


It  is  assumed  that  the  boundary  or  terminal  conditions  are 
specified 


o      o 


X 


Y„  <tJ  =  Y 
o     o  o 


o     o 


o 


x  (t  )  =  i 

o  o 

Y  (O =  Y^ 

o  o 

Z  (t   )  =  Z 
o  o 


X^)  -  xx 
Y1(V  =  Yl 

zx(tx)  =  zY 


Xj  (tx)  =  Xj 

Yi  <V  -  Yi 

Zj  (tx)  -  zx 


The  objective  is  to  select  the  functions  x(t),   y(t),    and  z(t)  so 

i       j,      i  Mo 

that  the  fuel  ratio-r-—    is  a  minimum  (i.e.,   we  desire  initial  fuel 

mass,    M     as  small  *  as  possible  and  still  attain  our  objective 

subject  to  the  boundary  conditions). 

The  minimization  of  the  above  integral  is  a  problem  in  the 
calculus  of  variations.  Its  solution  is  highly  detailed  and  is  not 
presented  here.     The  conclusions  and  interpretations  of  existing 
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solutions  are  helpful  and  interesting,   however,   and  are  stated 

37 
here  without  proof  . 


These  conclusions  are  due  to  D.  F,   Lawden,    one  of  the 
more  active  writers  in  the  field. 

Let  us  consider  a  rocket  orbiting  the  earth  circularly.     The 
optimal  path  is  desired  for  travel  from  this  orbit  to  entry  into 
circular  orbit  about  another  planet.     The  conclusions  of  Lawden 
indicate  that  an  optimal  path  can  only  be  followed  in  a  vacuum  by 
the  application  of  certain  impulsive  thrusts  from  the  rocket  motor, 
the  rocket  coasting  freely  under  gravity  during  the  intervals 
between  such  thrusts.     These  coast  periods  are  essentially  null- 
thrust  arcs     and  are  mihimal  in  the  sense  that  neighboring  arcs 
can  be  followed  only  by  an  actual  expenditure  of  fuel.     To  establish 
an  absolute  minimal  mode  of  transfer  for  a  definite  planet  con- 
figuration,   cases  for  various  numbers  of  impulses  must  be  cora- 

35 
puted     and  the  minimum  selected.     Lawden  suggests  that  two 

impulses  usually  lead  to  the  optimal  solution.     It  should  be  re- 
membered, however ,   that  there  will  be  a  large  number  of  optimal 
solutions  corresponding  to  different  planet  positions.     Only  one  of 
these  will  be  an  absolute  optimal  solution  for  all  planet  con- 
figurations ,      The  absolute  optimal  solution  for  an  Earth-Mars 

flight  has  been  known  theoretically  for  some  time  if  certain 

22 
fundamental  assumptions  are  made.     The  analysis       is  that  of 

Hohmann,   who  concluded  in  192  5  that  the  problem  of  transferring 

a  rocket  from  one  circular  orbit  into  another  coplanar  with  the 

first  and  about  the  same  attracting  body  would  be  most  economical 

if  effected  along  an  elliptical  trajectory  tangential  at  its  apses  to 

both  circular  orbits.     Fuel  would  be  expended  rapidly  at  each  apse 

to  carry  the  rocket  into  and  out  of  the  transfer  ellipse.     A  state  of 

free  flight  would  exist  between  the  two  impulse  trajectory.     This 

classic  development,   the  so-called  Hohmann  ellipse,   has  been 

verified  by  Lawden  and    is  consistent  with  the  remarks  above.     Note 

that  the  optimum  path  as  developed  and  computed  in  this  paper  on 


101 


the  basis  of  a  minimum  initial  velocity  is  in  fact  the  Hohmann 
ellipse.     The  importance,   however,    of  knowing  the  optimum  for 
other  positions  of  Mars  and  Earth  so  that  flight  can  be  initiated 
often,    is  not  to  be  underestimated. 


Mars 
Orbit 


Rocket  Path 


Decelerating 
Impulse 


Accelerating 
Impulse 


Fig,    6.4    The  Theoretical  Mode  of  Transfer 


6,  7    Rocket  Thrust  Direction 

Ordinarily,   the  mode  of  directing  the  thrust  which  results 
in  maximum  total  energy  is  found  to  be  alignment  of  the  effective 
exhaust  velocity  with  the  velocity  vector  of  the  rocket  (i.e.  align 
the  direction  of  thrust  with  the  tangent  to  the  trajectory.     This 
will  ensure  that  the  rate  of  increase  of  total  energy  is  maintained 
at  its  maximum  value  at  every  instant .     This  program  will  be  an 
economical  one  since  the  object  of  the  maneuver  is  to  raise  the 

total  energy  as     rapidly  as  possible  in  order  to  achieve  escape 

34 

and  sufficient  velocity      However,    it  has  been  proved       that  this 

program  may  not  be  an  absolute  optimal  method,   although  close 
to  it.    Under  some  circumstances,    when  departing  a  circular  orbit, 
it  is  not  economical  to  maintain  the  rate  of  increase  of  energy  as 
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its  highest  value  over  the  entire  trajectory  of  escape.     The 

author  claims  that  after  having  escaped,   it  is  advantageous  first 

to  direct  the  thrust  so  as  to  oppose  the  motion.     The  rocket  then 

drifts  toward  the  center  of  attraction  along  an  elliptical  arc.     At 

its  closest  point  of  approach,   an  impulsive  thrust  is  generated  in 

the  direction  of  motion.     The  velocity  is  now  large,   and  the  total 

energy  increases  at  a  rapid  rate.     The  rocket  is  transferred  into 

a  hyperbolic  orbit  along  which  it  coasts  to  its  destination  and  the 

overall  fuel  consumption  is  found  to  be  less  than  if  escape  had 

been  achieved  by  application  of  a  single  tangential  thrust.    The 

reader  is  referred  to  reference       for  information  as  to  how 

advantage  can  be  taken  of  such  maneuvers  for  general  conservation 

of  energy.     A  good  approximation  to  the  optimal  programming  of 

thrust  direction  under  departure  from  circular  orbit  may  be 

achieved  by  maintaining  alignment  between  directions  of  thrust 

and  motion, 

51 
Other  writers       consider  it  an  empirical  rule  that  a  journey 

from  one  to  another  stable  orbit  round  a  central  body  will  be 

completed  in  the  shortest  time  only  if  power  is  applied  in  two 

impulses,   one  at  beginning  and  one  at  end  of  journey. 


CHAPTER  7 
NAVIGATION  IN  SPACE 

A,    POSITION  IN  SPACE 

7. 1      General  Discussion  of  Problem 

Mid-course  interplanetary  navigation  differs   quite  markedly 

from  terrestrial  position  determination  due  primarily  to  the  need 

to  obtain  a  fix  within  a  volume  instead  of  upon  the  surface  of  a 

sphere  (Earth)  of  known  radius.     Terrestrial  positions  can  be 

obtained  by  simple  spherical  trigonometric  relations  and  can  be 

9 
further  simplified  to  a  problem  in  plane  geometry  and  time. 

Although  it  is  obvious  that  the  problem  of  Position  in  Space 
is  three  dimensional,    several  writers  have  postulated  that  a  suf- 
ficiently accurate  solution  can  be  obtained  by  considering  the 
planets  to  be  coplanar  and  working  with  only  two  dimensions  in 

O       CO 

the  plane  of  the  ecliptic.    '  Part  of  the  argument  for  this 

approach  is  found  in  star  tracker  limitations,    errors  in  astro- 
nomical measurements,   the  assumption  that  terminal  guidance 
will  correct  for  errors  incurred  during  mid-course  transit,    and 
the  relatively  small  angles  of  inclinations  of  the  planet  orbits 
from  the  ecliptic  plane. 

7.  2      Factors  Indicating  Necessity  for  "On  Board"  Automatic  Solution 

If  we  postulate  that  accurate  positional  data  is  needed  during 
mid -course  transit,    several  factors  indicated  that  the  problem 
must  be  solved  on  board  the  vehicle  and  by  automatic  computers 
of  some  sort. 

1.    Due  to  high  vehicle  velocities  (in  the  neighborhood  of 
1,  200  miles  per  min)  time  delays  of  any  sort  should  be 
avoided  wherever  possible.     This  eliminates  manual 
computation  and  sight  taking  in  the  vehicle  itself  or  on 
the  earth . 
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2.  Radio  frequency  transmission  in  space  has  a  rather 

o 

sizeable  time  delay,    (A  10     mile  round  trip  signal  requires 
approximately  18  minutes) 

3.  Radio  frequency  signal  interference  in  space  is  un- 
certain and  apparently  unpredictable 'when  sent  to  or 
from  Earth  due  to  the  Earth's  atmosphere  and  extraneous 
radio  frequency  waves  in  space  itself, 

4.  Power  requirements  for  dependable  radio  trans  - 
missions  during  mid-course  appear  to  be  far  too  great  for 
installation  on  any  vehicle  within  the  near  future . 

5.  Optical  tracking  of  the  vehicle  from  earth  is  not  possi- 

69 
ble  with  any  foreseeable  astro  telescopic  equipment. 

6.  Radar  tracking  of  the  vehicle  from  earth  is  ruled  out 
by  2  and  3 . 


7.  3     Selection  of  a  Reference  Coordinate  System 

Celestial  bodies  may  be  located  with  respect  to  several 
different  reference  systems  which  are  shown  in  Fig.  7.1       and 
Table  7.1  ,     For  mid-course  interplanetary  position  de- 

termination two  of  these  systems  can  be  discarded  immediately. 
These  are  the  Geocentric  Horizon  and  Equatorial.     Coordinate 
systems  using  Galactic  quantities  are  not  particularly  convenient 
for  travel  limited  to  our  solar  system.     Of  the  remaining  only  the 
heliocentric  ecliptic  system  will  be  considered  for  use  in  this 
paper.     This  system  has  the  following  features  which  lead  to 
this  choice. 

1.  The  reference  zero  latitude  plane  includes  the  Earth 
and  Sun  at  all  times. 

2.  All  other  planet's  orbits  lie  close  to  this  plane, 

3.  NumerousEphomerides  are  available  based  in  this 
system  simplifying  programming  of  any  memory  drum  for 
a  computer 
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7 .  4      Establishing  the  Reference  Coordinate  System 

It  appears  that  the  reference  system  orientation  will  have  to 

be  obtained  by  use  of  the  fixed  stars,  an  extremely  sensitive  gyro 

71 
inertial  platform,       or  a  combination  of  the  two.     The  use  of 

gyroscopic  elements  by  themselves  for  establishing  and  main- 
taining a  chosen  reference  coordinate  system  to  a  sufficient  degree 
of  accuracy  to  determine  accurate  angular  measurements  and  rates 
in  space  appears  to  be  beyond  engineering  capabilities  at  present. 
All  instrument  gyros  have  drift  uncertainties  to  some  extent  and 
very  definite  lower  limits  of  sensitivities.     To  be  useful  for  ac- 
curate interplanetary  guidance  the  threshold  of  sensitivity  would 
have  to  be  in  the  order  of  10"   g.     Since  it  is  necessary  to  con- 
template transit  times  of  over  100  days,   it  appears  that  star 
trackers  locked  on  the  fixed  stars  will  have  to  be  used  to  a  cer- 
tain degree  in  either  assisting  a  gyro  system  or  replacing  it 
entirely.     At  certain  periods  of  the  transit  it  may  be  advantageous 
to  use  gyroscopic  elements  as  a  memory  circuit  or  coast  device 
while  a  star  is  temporarily  obscured  or  its  line  of  sight  becomes 
nearly  coincident  to  that  of  a  stronger  source  of  light  such  as  a 
planet  or  the  sun.     It  also  appears  that  during  periods  of  radical 
maneuvering  or  high  thrust  a  gyroscopic  device  may  be  needed 
to  insure  holding  track  on  a  desired  star. 

If  we  postulate  the  use  of  photoelectric  star  trackers,    es- 
tablishing the  desired  reference  coordinate  system  can  be  ac- 
complished by    using  only  two  fixed  stars  approximately  90     apart 
in  the  same  meridian.     However,   three  stars  near  quadrature  to 
each  other,   or  three  stars  near  the  chosen  zero  latitude  spaced 
approximately  120     apart,   will  establish  a  more  stable  and  secure 
reference.     Monitoring  three  fixed  stars  would  allow  for  the 
possibility  of  one  of  them  becoming  momentarily  obscured  or 
temporarily  loosing  tracking  signal.     Section  8,  64  discusses 
the  accuracy  requirements  for  such  systems. 


109 


7,  5     Determining  Celestial  Latitude  and  Longitude  of  the  Vehicle 

After  establishing  the  reference  coordinate  plane  and  meridian, 
the  celestial  latitude  and  longitude  of  the  vehicle  can  be  determined 
from  the  line  of  sight  to  the  sun  as  seen  from  the  vehicle,     This 
is  illustrated  in  Fig.  7.2  and  as  shown,   the  vehicle  heliocentric 
latitude  is  equal  to  the  negative  of  the  angle  between  the  reference 
zero  latitude  plane  and  the  vehicle  sun  line  measured  in  a  meridian 
of  the  reference  coordinate  system.     The  heliocentric  longitude  of 
the  vehicle  is  equal  to  the  angle  between  the  reference  zero  longitude 
and  the  meridian  containing  the  vehicle  sun  line  measured  in  the 
direction  of  rotation  of  the  planets,   plus  180    . 

7.  6     Determining  the  Radius  Vector  to  the  Sun 

After  establishing  the  basic  reference,  and  determining  the 
vehicle  latitude  and  longitude  with  respect  to  the  sun,    it  is   then 
necessary  to  obtain  the  distance  of  the  vehicle  from  the  sun.     This 
can  be  accomplished  using  only  one  planet  tracked  with  respect  to 

the  established  reference  system  by  a  simple  method  proposed  by 

w     *        +   67 
Vertregt, 

We  can  obtain  from  a  memory  drum  or  an  ephemeris  of 
heliocentric  ecliptic  positions  of  the  selected  planet  as  a  function 
of  time; 

1  Longitude  of  the  planet  {  X       ) 

2  .    Latitude  of  the  planet      (  L      ) 

v  p  v 

3,    Distance  of  planet  from  sun  (   R       ) 

To  obtain  the  radial  distance  of  the  vehicle  from  the  sun 
referring  to  Fig.  7.3,   we  can  proceed  as  follows: 
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Fig.    7.2     Determination  of  heliocentric  latitude  and  longitude 
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Fig.    7.3     Geometrical  relationships  determining  radial 
distance  from  *\\c  sun 
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If  the  earth  is  one  of  the  planets  used  for  determining  positon, 
and  we  use  the  plane  of  the  ecliptic  as  the  zero  latitude  reference 
plane,   the  solution  is  somewhat  simplified  due  to  the  fact  that 
earth  heliocentric  ecliptic  latitude  is  at  all  times  zero.   (Fig.  7.  4  ) 
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0)AA 


Fig.    7.4     Geometrical  relationships  determining  radial  distance 
from  the  sun  with  earth  as  the  reference  planet 


114 


The  determination  of  range  of  the  vehicle  from  the  sun  can 
be  made  using  only  one  planet  after  L.  and   X.  are  known. 
However,    since  we  know  appreciable  errors  exist  in  the  values 
given  for  planetary  positions  as  shown  in  Section  8.  62,   the  use 
of  two  or  more  planets  should  tend  to  narrow  the  error  down. 

7.  7      Some  Additional  Methods  for  Position  Determination 

A  second  method  for  angular  determination  of  the  sun  and 

fi  7 

planets  is  photographic  u>  /).     .  By  the  taking  of  simultaneous 

pictures  of  a  given  small  area  (such  as  5   )  centered  about  the 
selected  body,   the  position  of  the  body  with  respect  to  the 
background  of  recognizable  fixed  stars  can  be  obtained  by  taking 
measurements  from  the  developed  plate.     This  does  not  appear 
practical  even  for  a  manned  vehicle,    due  to  time  lag,   the  need 
for  continuous  or  very  frequent  data  to  effect  smooth  system 
performance  and,   in  the  case  of  a  manned  vehicle,   an  extra 
burden  upon  the  crew. 

Similarly  it  has  been  suggested  that  the  angular  measurement 
of  the  diameter  of  solar  bodies  might  be  used  for  range  determi- 
nation.    This  could  be  photoelectric  or  photographic  in  nature. 
Both  these  systems  do  not  appear  practical  except  possibly  for   . 
terminal  guidance  where  the  diameter  of  the  body  being  approached 
becomes  appreciable. 

Further  postulated  ranging  methods  are  (1)  measurement  of 

49 
incident  solar  radiation  intensity,        (2)  measurement  of  reflected 

solar  radiation  from  planets,    and  (3)  radar  ranging. 

Some  basic  arguments  against  use  of  the  above  methods  for 
mid-course  ranging  are  as  follows: 

1.    The  variation  in  intensity  of  solar  radiation  does  not 
seem  to  be  great  enough  for  any  but  the  most  coarse  range 
determination      Using  Fig.   2.  2        we  see  that  for  a  trip 
between  Earth  and  Mars  this  variation  in  radiation  is 
approximately  1  erg  per  square  inch  per  second  for  a 
change  in  range  of  one  mile. 

115 


2.  Reflected  solar  radiation  also  has  the  same  general 
gradiant  as  above  in  (1)  plus  the  uncertain  and  possible 
widely  varmg  amounts  that  might  be  reflected  from  a 
body  do  to  changes  in  planetary  atmospheric  condition 
and  orientation  of  lighted  surface  relative  to  the  vehicle. 

3.  Radar  ranging  has  the  same  disadvantage  as  listed  in 
Section  7.  2 

7.  8      Determination  of  the  Shape  of  Actual  Transfer  Ellipse 

If  we  assume  that  after  obtaining  a  series  of  acceptable 
vehicle  positions  we  want  to  know  the  elements  and  orientation  of 
the  ellipse  being  traversed,    we  can  proceed  after  first  converting 
the  differences  in  latitude  and  longitude  of  three  selected  positions 
into  angular  displacements  along  their  great  circle  locus  (6^  -  6,). 
This  can  be  accomplished  by  standard  coordinate  transformation. 

Determination  of  three  of  the  basic  elements  of  the  unpowered 
orbit,   6  ,    e,    and  a,    can  be  accomplished  by  the  use  of  three  se- 
lected positions  of  the  vehicle  with  respect  to  the  sun.     For  any 
reasonable  degree  of  accuracy  the  time  interval  between  selected 
positions  may  have  to  be  several  days  .     A  method  of  mathematical 
reduction  which  can  be  easily  programmed  for  a  computer  has  been 
suggested  by  Vertregt       and  is  adapted  herein  to  the  nomenclature 
of  this  paper  in  Derivation  Summary  7.1. 

Two  assumptions  are  made  in  this  method  which  are  common 
to  many  proposed  non-powered  interplanetary  mid-course  tra- 
jectories     These  are. 

1.    Vehicle  is  following  an  unperturbed  Keplerian  ellipse, 

2  .    Sun  is  only  gravitational  force  considered,    and  is  one 
of  the  foci  of  this  ellipse. 

From  Derivation  Summary  7,1  we  obtain 
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DETERMINATION  OF  THREE  BASIC  ELEMENTS  OF  ACTUAL 
ELLIPTICAL  TRAJECTORY  FROM  OBSERVED  DATA 


The  following  method  of  determining  0,   e,   and  a  is  adapted 
from  reference  67  to  the  nomenclature  of  this  paper. 

Using  a  basic  equation  of  any  ellipse  and  referring  to  Fig.  7.  5 
we  obtain  the  following: 


a  (1  -  e2) 

RA    = 

i 
1+  e  cos  6, 

a  (1  -  e2) 

RA       = 

A2 

i 
1+e  cos  #2 

a(  1-  e2) 

RA      = 

A3 

i 
1  +  e  cos  Or. 

(5) 
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(6) 


(7) 


Combining  (5)  and  (6)  and  (5)  and  (7)  we   get 

i  i 

RA    +  RA    e  cos  6,  =  RA     +  RA     e  cos  6n  (8) 

Al         Al  l         A2         A2  2 

RA   +  RA    e  cos  6,  =  RA    +  RA      e  cos  el  (9) 

Al         Al  X        A3         A3  3 

Rearranging  (8)  and  (9) 

Ra     -  Ra    =    e  (Ra     cos  0i  "  Ra     cos  6n  ) 

A2         A1  A1  1         A2  2 

RA9  "  RAt 
e= -t * r  (10) 

RA     cos  6,  -  RA      cos  6n 
Al  l         A2  2 
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R 


R 


=  e  (RA     cos  0,  -  RA     cos  0Q  ) 
Al  Al  l         A3  3 


R  A  "      R  A 

A3         Al 


^  A  "     ^A 

A2  Al 


RA    cos0,~RA     cosl      RA    cos0,-RA    cos  0o 
Ax  1       A3  3        Ax  1       A2 


(11) 


From  (11)  we  obtain 

t  i 

RA  (RA     -RA    )  cos  0.  -  RA     (RA    -  RA    )  cos  0O 

Al     A2        A2               X  A2       Al         A3               2 


■b  RA      (RA     -  RA      )  cos  0o 
A3       Al         A2  3 


\-  V(V  V 


v<vS> 


0 


(12) 

(13) 
(14) 


(0„  -  0,  )  and  (0„  -  0,  )  have  been  determined  along  with  R.     , 

R»      ,   R.      by  establishing  the  three  instantaneous  positions  of  the 
vehicle,     we  can  now  substitute  (13)  and  (14)  into  (12) 


RA    (RA     -  RA    )  cos  0.  -  RA     (RA  -  RA    )  cos 

Al       A2         A3  l         A2       Al       A3 

r  i 

, +  RA     (RA    -  RA     )  cos      0,  +  (0O  -  0,  ) 
A3       Al         A2  L  l         3         l  - 


e[He2-e[) 


=    0 


(15) 


Expanding  the  summed  angle  terms  in  (15)  leads  to 

RA    (RA     -RA    )  cos  0,  -  RA    (R.    -RA    )   cos  0.  cos  (0O-0J- 
Al      A2        A3  1         A2      Al        A3  L 


i  i         i 

-  sin  0,  sin  (0    -  0  ) 

i  it 

-  sin  0,  sin  (0„  -  0   ) 


+  RA    (R      -RA    ) 

A3      Al        A2 

=   0 


i  i         t 

cos  0,  cos  (0o-  0,  ) 


(16) 
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Solving  (16)  for  0,  we  get 


6    «    tan 


JRA1(RA2  "  RA3>  "  RA2(RA1-  RA3>  COS  (e2  "  flI>  +  W  %>  c°8  <63  "  9l* 


-R 


A2(RA1-  RA3>  Sin  (62  "  61)+RA3(RA1-  RA2>  sin  <63  "  V 


(17) 


From  (11),    substituting  the  value  obtained  in  (17)  we  obtain 


e    =■ 


A3         Al 


RA     cos  0,  -  RA     cos  6n 
Al  l         A2  3 


(18) 


Substituting  the  value  obtained  in  (18)  into  (7)  and  rearranging 
leads  to 


a  = 


R.      (1+  e  cos  0  ) 

:     2 

1  -  e 


(19) 
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=    tan 


R.    (RA      -  R.    )  -  R.    (RA   -  R.    )  cos  (6„  -  6,)+  RA    (R  .    -  RA    )  cos  (6 


-i       *i    A: 


2  "*3 


A2      *1 


A3  ^1        V 


"A2(RA1- 


RA    )  sin  (i 
3 


ei)+RA3(RA1-RA2,Sm<e3"e,) 


R 


R 


e  = 


vl 


g  i 

RA     cos  0,  -  RA     cos  0O 

Ax  1         A2  3 


R  ,      (1     e  cos  0Q 
^3 1 

i         2 
1  -  e 


(7=5) 


(7-6) 


(7e-)7) 


7.  9      Determination  of  Orientation  of  Actual  Transfer  Ellipse 

Two  additional  elements  of  the  transfer  ellipse,  fi  and  i,   are 

required  for  determining  the  orientation  of  the  vehicle  trajectory 

with  respect  to  the  reference  plane.     These  can  be  obtained  by 

fi  7 
applying  the  method  suggested  in  reference        using  two  positions 

of  the  vehicle.     A.  rather  obvious  observation  can  be  made  that  the 

further  these  two  positions  are  from  each  other  the  more  accurate 

will  be  the  values  of  Q  and  i  obtained.     However,    since  corrective 

action  is  more  economical  if  performed  early  in  the  flight,    it  is 

apparent  that  Q,  and  i  should  be  determined  as  early  as  possible. 

If  the  trajectory  is  satisfactory  at  this  time  more  accurate  values 

can  be  obtained  at  the  next  determination. 

From  the  first,  vehicle  position,   as  illustrated  in  Fig.  7.  6    we 
can  obtain  the  following: 

to  NN* 

,o 


Vs 


A,. 


R 


1 


1  '© 


Z.  -   RA 

i  A 


sin  L 


1 


90 


1 


n 


i 


R 


cos  A. 


RA^S1}0     RA}  COS 


90    -  ij^R^    sin  i^ 
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therefore 


sin  L 


sin    1    = 


Ai 


sin    r] 


(7=8) 


For  the  second  position  we  can  likewise  determine  that 


sin  L 


sin    1  = 


sin[n1+(e2  -6.) 


Combining  the  above 


sin  L 


sin    n 


sin  L 


sin  L 


tl   -   . .  

1  sm[n1  +  (0,'-  6')]       sin  n  cos  {  02  -  0')+  cos  n  sin(02~  &\) 


Multiplying    by  sin  n,  and  regrouping  leads  to 


r?1  =  tan 


sin  L  .     sin  idn  -  0,) 
A  z         i 


sin  L.      -  sin  L,     cos  (0„  -  ©Jj 


(7=9) 


n,  can  now  be  converted  to  the  plane  of  the  ecliptic.     If  we  then 
subtract  it  from  A.,    we  have    ft. 

After  determining    rj,  we  can  then  solve  for  the  inclination  of 
the  orbit  of  the  vehicle  using  equation  (7=8) 


(7-10) 


i  ; 

.     -1 

=  sin 

"sin  L  A  " 
A! 
sin  7], 
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B,   IN-FLIGHT  CORRECTIONS 


7,10     Introduction 

It  is  considered  inevitable  that  the  rocket  may  fail  to  follow 
its  prescribed  path  in  space      Errors  in  initial  conditions  or 
unknown  gravitational  influences  will  cause  deviations  from  the 
predicted  trajectory. 

As  soon  as  differences  between  actual  position(as  described 
in  Part  A  above  and  predicted  position  become  known, 
a  new  situation  arises.     The  decision  must  be  made  to  either 
return  to  the  previously  calculated  orbit  by  means  of  corrective 
thrust,    or  to  recompute  a  new  path  which  will  be  within  the 
capabilities  of  the  rocket.     Time  is  an  important  consideration  in 
this  discussion  since  corrective  action  must  place  the  rocket  on 
a  path  which  will  result  in  almost  simultaneous  coincidence  of 
Mars  and  rocket.     This  is  discussed  further  below. 

It  is  presumed  that  almost  continuous  position  finding  data 
will  be  available  so  that  comparisons  can,    in  fact,   be  made 
between  actual  and  predicted  position.     It  is  also  reasonable  to 
assume  that  deviations,   hence  corrections,   will  be  incremental 
in  nature.     For  the  sake  of  simplicity,    let  us  also  assume  that 
computation  of  a  new  trajectory,    when  required,    can  take  place 
instantaneously.     Under  actual  conditions,    solution  time  can  be 
allowed  for  with  no  great  difficulty. 

The  desirability  of  placing  the  rocket  on  a  new  path  rather 
than  forcing  it  to  adhere  to  a  previously  defined  path  has  been 
mentioned.     The  latter  method  implies,    in  effect,   the  application 
of  continuous  corrective  thrust.     Reference  26  has  indicated  an 
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analytical  approach  to  the  thrust  requirements  necessary  to  fly 
a  prescribed  path.     Until  such  time  as  fuel  supply  poses  no 
problem,   however,    continuous  corrective  thrust  is  prohibitive. 
As  little  application  of  corrective  thrust  as  possible  is  desired, 
since  there  is  little  available  for  expenditure. 


» 


The  objective  of  corrective  action  is  not  necessarily  to 
achieve  a  pinpoint  interception  of  Mars,   but  rather  to  reach  a  posi- 
tion favorable  for  entry  into  orbit  about  Mars.     This  might 
conceivably  be  a  position  at  some  distance  thousands  of  miles 
from  the  planet  as  well  as  on  the  planet  itself,    since  the 
gravitational  influence  of  the  arrival  planet  and  decelerating 
thrust  will  combine  to  make  entry  into  orbit  a  less  significant 
consideration.     Therefore,    some  error  can  be  tolerated.     However, 
when  and  if  deviations  occur  which  are  outside  tolerance  limit, 
corrective  thrust  becomes  necessary.     Obviously,   the  tolerable 
deviation  is  smaller  at  the  early  stages  of  flight. 

For  purposes  of  analysis  let  us  assume  that  the  allowable 

tolerance  has  been  exceeded  and  that  corrective  action  is  called 

for.     The  rocket  is  moving  primarily  under  the  attraction  of  a 

single  body,   the  sun,   which  is  justifiable  in  accordance  with  the 

calculations  of  section  3.6,     We  will  also  consider  Mars  as  a 

massless  point  to  be  "intercepted"  (in  the  sense  defined  above) 

at  a  specified  time. 

28 
The  subsequent  analysis  parallels  that  of  Lawden 

although  different  in  many  respects. 

7.11     Analytical  Development 

i 
Let  r     ,  0     define  present  position  in  space  at  time  t     .    Let 

r,  ,    0,  define  the  position  the  rocket  would  occupy  at  time  t     had 

flown  the  correct  path.     r„  ,    69  define  the  position  of  Mars  at  time 
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t2,  the  planned  time  of  interception  It  is  desired  to  intercept 
Mars  at  time  to  via  a  route  which  obviates  the  necessity  of  re- 
turning to  the  original  trajectory.  Figure  7,  7  is  helpful  in  the 
interpretation  to  follow,  The  route  is  a  conic  passing  through 
P  and  J. 


Mars 
Trajectory 


/ 


to 

L 


o 


r:^e 


,ic 


x0. 


reference  line 


perihelion  of  predicted 
path 


Fig.    7,  7     Interception  at  Originally  Predicted 

Time,     t« 


Clearly,    it  is  a  restriction  to  be  forced  to  intercept  the 
planet  at  time.    t?,   but  results  in  a  great  simplification  in  the 
computer  requirements.     Let  us  pursue  this  analysis  first. 

The  conic  path  to  be  flown  is  presumed  to  be  a  segment  of 
an  ellipse  (because  of  velocity  requirements)  with  sun  as  focus 
The  line  of  apsides  of  this  new  ellipse,    whose  elements  are 
designated  (a«,    e?)  is  aligned  at  an  angle   u  with  the  reference 
line  of  the  originally  prescribed  ellipse  (elements  a,,e,) 

From  the  fundamental  elliptical  equation  4-1 


a2  (1  ~  e2   ) 


1  +  e0  cos  (6      -   u  ) 
2.  o 


(7-11) 
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a„  (1  -  e0    ) 

r2    = i—^—  (7-12) 

1  +  e2  (0„  -  to  ) 


where  r   ,    r9,  6    ,    and  90  are  known. 

O  £         o  ^ 


Position  in  an  elliptical  orbit  at  any  time  t  can  be  found  from 

47 
the  relationship 


? 


2     .         '  i  e,  +  cos  0, 


2?rt      eVl-e     sin    0i       D;«-l_l L^Z 

2 


1  -  sin      f  +  -s- 


i  +  ei  cos  e[  l  +  e  cos  6\ 


where  t  is  measured  from  t  =  0  which  occurs  at  perihelion  passage 

of  ellipse  being  considered. 

44 
T,   the  period  of  the  ellipse^can  be  expressed       as  a  function 

of  a,    and  is  usually  known, 

2-23/2 
T2  '   ^UW  (7'13) 

From  these  relationships, 


[~         2  '  ' 

27r  t2        eo'^1  "  e2       sin  ^  ®2  "   u^  -1  e2  +  cos^  2  ~  u  '        7r 

rp  j  Sill  I  1      q" 

2  1  +  e2  cos  (  02  -  to  )  l+e?  cos(02  -  co  ) 

and  I g  >  » 

27T  t         e0v  1  -  en     sin  ( 9      » to)  -,  e0+  cos  (6    -  to) 

O  £  6  O  -1       6  O  7T 

=  j -    sm      j — -    +y 

2  1  +  eQ  cos  (0     -to  )  1+  en  cos  (0O  -to) 

where  tn  and  t     are  unknowns,    measured  from  the  reference  line 

<£  O 

inclined  at  an  unknown  angle   to    to  a  known  reference  line . 

If  the  time  of  travel  from  r     to  rn  is  constrained  to  be  the 

o  I 

originally  planned  interval  of  time  from  r,  to  r?,   the  interval 

t2  -  t     becomes  a  known  quantity,    and  t?  and  t     need  not  be  known 

individually. 


12« 


T 

tn     ~    t        =     At    =    7T- 

2        o  2tt 


I  2  '  /         2       ' 

e2  2      sin(02-u)      e2\/l-e2    (6Q-u) 


i  i 

1  +  erj  cos  ( 0„  -    w)  1  +  eQ  cos  (0    -to) 

2  2  ^  o 


i  i 

=1  e2  +  cos  (02  -co)  i    e2+cos(0Q-  to) 

-  sin  j -i-sin      , 

1+  e2  cos  (02  -  to)  1+  e2cos(0    -to) 

(7-14) 


In  summary,    position  finding  data  will  provide  r    ,  0    .     The 
quantities  r,  ,  6.    and  r„,    Or,  are  known  in  advance  as  functions  of 
time.     Thus,  At     the  time  till  interception,    is  known  for  any  r, , 
6y     Equations  7   11     7-14involve  four  unknowns  a?     e2,  to,    and  T2, 
and  are  solvable  by  iterative  methods 

Since  corrections  are  assumed  to  be  incremental  in  nature 
when  compared  with  the  elements  of  the  path  being  flown(a   ,  e    , 
to  ),    a?,    e~  and  to,   may  be  replaced  by       a  +  6a,    e  +  6e,  and 
to  +    6to.     It  is  assumed  T?  has  been  replaced  by  its  value  in 
terms  of  a? 

If  these  substitutions  are  made  and  small  angle  assumptions 

are  valid  (sin  6to   -   6u>,    cos  6to  =  1)  equations  in  the  three  unknown 

6a     6e,   and  6to  result      a    .    e     and  to     (assumed  zero,    initially) 

o       o  o  J 

can  be  found  as  described  in  the  navigational  section  and  are  good 
first  approximations  to  a      e    and  to  . 

Successive  improvements  of  the  accuracy  of  6a,    6e,    and  6to 
can  be  made  if  desired 

7„12   Corrective  Action 

Since  changes  must  inevitably  be  reflected  in  an  alteration  of 
velocity  and  angle  of  fire,    we  may  again  write  two  basic  re- 
lationships 

(v     +  Av)2  =GM(JL    _i_) 
o  r2        a2 
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! 

e9  sin  (  6  9  -  to9) 

tan  ( </>    +  A<f>     ) 


1  »  e9  cos  (02  -  402 ) 

Replacing  a«,    e9  and  to-  by  a+  6a,    e9  +  6e  and  u+  6u,   we 
arrive  at 

(v    +Av)2     ,    GM   [    |- 


o  r9        a  +  6a 


Av  = 


GM  ' 


r9       a  +  6a 


2 

„v 
o 


where  v     is  the    velocity  at  the  time  corrective  action  is  deemed 

o  J 

necessary,    and  <j)     is  the  corresponding  angle  of  orientation  as 

defined  earlier.     Similarly,    direction  must  be  altered  to  effect 

Ad     which  can  be  written 
yo 

i 
,    (e  -t-6e)  sin  (09  =  u   -  6w) 
A0     =  tan j -  <j) 


o 


l-(e+6e)  cos  (02=w"^w^ 


a  +6a,    e  +  6e,    and  u    +  6u  are  presumed  to  be  known  accurately 
as  described  in  the  preceding  paragraphs. 

Corrective  thrust  must  be  actuated  to  produce  the  changes 
required. 

7.13    Interception  at  Newly  Calculated  Time 

Let  us  now  consider  the  prospect  of  intercepting  Mars  at  some 
time  other  than  the  originally  planned  time. 
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Mars 
Path 


elements 
&2,  e2 


reference  line 
i  predicted 


Fig     7,8    Interception  at  Newly  Calculated 

Time 


It  is  recognized  at  t     that  present  position,    r    ,  9[    is  un- 

&  0  r-  f  O        O 

satisfactory,  i.  e     a  continuation  along  present  track  will  not 
achieve  interception.     Moreover,    assume  that  rocket  corrective 
thrust  capabilities  will  not  permit  an  interception     at  r?J    0L.   the 
originally  planned  arrival  point.     We  desire  a  conic  (elliptical) 
path  from  r       6[    to  some  arbitrarily  selected  Mars  position,   r^, 
6>'o  which  has  been  chosen  because  it  is  within  the  rocket's  thrust 
capabilities      This  new  position  may  also,    in  fact,    represent  the 
position  of  Mars  which  can  be  most  economically  reached.     The 
choice  of  this  position  does  present  a  formidable  problem,    since 
the  optimization  process  is  necessarily  involved.     The  selection 
of  a  Mars  position  is  not  discussed  here  although  its  importance 
is  emphasized. 
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Presumably,   there  is  on  board  a  means  of  establishing  Mars 

position  as  a  function  of  some  base  time  -  a  Mars  Ephemeris  of 

sorts.     Fixing  the  desired  interception  point  will  fix  the  time,  t.„, 

of  interception,     t     is  known  from  an  accurate  time  piece.     Thus, 

the  interval  of     corrective  travel  becomes  *t0  -  t     ,     These  times 

3        o 

are  both  based  upon  a  common  reference.     The  time  interval  is 
therefore  determinable. 


Employing  the  previously  developed  reasoning  procedures, 

2, 


r     - 
o 


a   (1  -  e2 


l  +  en  cos  (6    +u) 
2  o 


r3" 


a2(l  -  ea*> 


l  +  e2  cos  (0„+  to) 


At  =  t. 


T 
2 

*o  '  2tt 


f- 


T 


e„     sm 


>+  jj 


1+  e2  cos  (6„  +w) 


/ 21  ' 

22Y1"e2     sin  ^Go+  u^ 

i 

1+  e9  cos  (6     +  u)  ) 


1                                                                 I 
i        e2+  cos  (6„  +  u)              ,      e2+  cos  (0     +  to  ) 
-  sin   '       ; fsin   '      i 


1+  e2  cos  {dr,  +u 


1+  en  cos  (6    +w  ) 
2  o 


and 


2w  a 


3/2 


y/^M 


By  the  methods  previously  described,   the  elements  of     the 
new  path  a„  ,   e~  ,  co9   can  be  found  and  corresponding  6a,    6e,   and 
6co's  determined.     It  is  now  a  relatively  simple  matter  to 
calculate  the  velocity  and  angular  changes  necessary  to  fly  the 
corrected  path. 


132 


7.14    Summary 


It  is  seen  that  both  approaches  are  essentially  the  same. 
From  the  standpoint  of  instrumentation,   however,    interception 
at  time  other  than  planned  requires  on  board  a  knowledge  of  Mars 
position  at  all  times  and  necessitates  the  selection  of  a  new  inter- 
ception point.     The  choice  of  this  point  is  essentially  an 
optimization  process. 


CHAPTER  8 
GUIDANCE  AND  CONTROL 


8. 1    General 

The  instrumentation  of  any  guidance  and  control  system  for 
interplanetary  travel  is  extremely  complex.     As  stated  previ- 
ously,  this  paper  studies  the  so-called  "combination"  trajectory 
(i.e.    -an  initial  thrust  period  followed  by  unpowered  free-fall 
flight;  modified,    if  neaessary,   by  corrective  thrusts).     Primary 
emphasis  has  been  placed  upon  what  has  been  defined  as  the  "mid- 
course"  of  this  type  trajectory.     Although  considerable  thought 
has  been  devoted  to  the  use  of  departure  and  arrival  orbits,    the 
applicability  of  this  paper  is  not  limited  solely  to  interorbital 
trajectories  .     Whether  the  vehicle  orbits  prior  to  commencing 
mid-course  travel  or  simply  continues  directly  from  the  ascent 
from  the  surface  of  the  Earth  does  not  affect  requirements  of  the 
vehicle  velocity  vector  orientation  and  magnitude  and  other  para- 
meters for  the  actual  transfer  period. 

This  Chapter  outlines  in  simple  form  a  suggested  Mid- 
course  Guidance  and  Control  system  based  upon  the  material 
developed  in  previous  portions  of  this  paper      The  major  problem 
of  reducing  the  system  to  the  least  possible  weight  and  size  has 
been  ignored  except  to  the  extent  that  simplification  has  been 
attempted  where  it  will  not  detrimentally  affect  the  overall  ac- 
complishment of  the  mission  (i.e.  -  the  arrival  of  the  rocket  at 
Mars  within  certain  limits). 

Time-sharing  of  a.  miniaturized  low-power  digital  computer, 
a  common  ephemeris  memory,    and  one  basic  time  generating 
mechanism  have  been  assumed.     This  appears  justifiable  since 
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the  elasped  time  between  computations  need  not  be  of  the  order  of 
fractions  of  a  second,   but  can  be  several  minutes  or  possibly  hours 
Digital  computation  is  necessary  due  to  the  extremely  low  rates  of 
change  of  rocket  position,   velocity,   and  acceleration;  and  the  need 
for  a  high  degree  of  accuracy. 

The  problems  of  instrumentation  have  been  divided  into 
several  sections  as  follows: 

(1)  The  Basic  System 

(2)  The  Tracking  Section 

(3)  The  Computing  Section 

(a)  Tracking  Line  Computation 

(b)  Prediction  Computation 

(c)  Comparing  Function 

(d)  Correction  Computation 

(4)  Thrust  and  Stabilization  Section 

(5)  Errors  and  Effects 

These  sections  are  treated  below 


8 ,  2  The  Basic  System 

The  functional  responsibilities  of  Guidance  and  Control  have 
been  segmented  into  three  basic  categories,   as  shown  in  Fig,    8-1. 
A  brief  general  description  follows 


Vhe  Tracking  Section 


The  Computing  Section 


Maintains  Lines -of -Sight  upon 
prescribed  bodies  and  provides 
Tracking  Line  data  therefrom  to 
the  Computing  Section. 

Performs  the  necessary  computations 
to  the  Tracking  Lines  to  obtain 
signals  proportional  to  the  cor- 
rections required  to  the  rocket 
velocity  vector  to  insure  ac- 
complishment of  the  mission. 
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The  Thrust  and  Stabi- 
lization Section  -    Provides  vehicle  orientation  and 

corrective  thrust  in  accordance 
with  the  received  signals , 

The  aforementioned  system  sections  will  be  examined  in 
more  detail  in  the  following  sections  of  this  paper.     The  Computing 
Section  of  the  rocket  involves  such  a  complex  problem  that  it  has 
been  subdivided  into  four  parts,    each  of  which  will  be  discussed  in 
its  own  Section  of  this  Chapter. 

8.  3  The  Tracking  Section 

The  tracking  section  envisioned  determines  tracking  lines 
from  the  vehicle  to  the  Sun  and  planets  measured  in  an  helio- 
centric ecliptic  coordinate  system.     These  angular  measurements 
are  then  sent  as  electrical  signals  (TL.)  to  the  tracking  line  com- 
puter in  the  tracking  section.    Two  systems  are  presented  which 
differ  in  method  of  space  reference  frame  determination,   number 
of  star  and  planet  trackers  required,    and  in  probable  degrees  of 
accuracy  and  reliability. 

8.  31    System  One 

The  first  system  obtains  space  reference  from  star  trackers 
rigidly  mounted  to  a  control  member  platform  isolated  from  vehi- 
cular motion  by  gimbals.     This  platform  is  "fixed"  in  space  and  is 
maintained  in  this  attitude  by  tracking  error  signals  generated  in 
the  fixed  star  tracking  units,   which  are  used  to  actuate  the  gimbal 
drive  motors.     A  gyro  package  is  mounted  on  the  controlled  member 
solely  to  assist  in  holding  track  during  periods  of  radical  maneuver 
and  to  aid  in  initial  lock  on.     Three  additional  star  trackers  are 
mounted  on  the  controlled  member,    each  with  its  own  tracking 
member  drive  system,   as  illustrated  in  Fig,   8.2  and  8.3.     One 
tracking  head  determines  the  rocket-sun  tracking  line  while  the 
others  determine  two  different  rocket -planet  tracking  lines. 
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SOLAR  SYSTEM 
TRACKERS 


-x- 


-^ 


MIDDLE  GIMBAL 
DRIVE  MOTOR 


SPACE  REFERENCE 
FIXED  STAR  TRACKERS 


OUTER  GIMBAL 
DRIVE   MOTOR 


TRACKER  VERTICAL  PLANE   >*,— 

DRIVE  MOTORS 


REFERENCE  PLANE 
CONTROLLED  MEMBER 


GYROS 


Fig.   8.2      Illustrative  schematic  of  tracking  system 

with  gyro  assist 
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Angular  measurements  of  these  lines  are  made  in  heliocentric 
ecliptic  coordinates,    (see  Fig,    7.1     )  and  transmitted  to  two 
separate  range  computers  within  the  tracking  line  computer  of 
the  computing  section 

8.  32  System  Two 

The  second  tracking  system  (Figs.  8,4  and  8.5)  consists  of 
a  space  reference  gyro  package,    and  only  three  star  trackers, 
The  reduction  in  the  number  of  trackers  used  is  made  solely  to 
simplify  construction  at  the  expense  of  system  accuracy      Three 
orthogonally  positioned  single  degree  of  freedom  gyros  are  as- 
sisted by  a  single  "fixed  star"  tracker  rigidly  mounted  to  the 
controlled  member      Smoothed  error  signals  from  this  tracker 
are  combined  with  the  signals  from  the  gyros,     This  resultant 
signal  is  then  used  to  actuate  the  gimbal  drive  motors  maintaining 
the  controlled  member  reference  plane  in  proper  orientation. 
Two  star  trackers  are  then  mounted  on  the  controlled  member, 
each  with  its  own  tracking  member  drive  system.     As  in  the  first 
system     one  tracker  determines  the  rocket-sun  line,    while  the 
second  tracker  determines  a  rocket -planet  line.     Angular  measure 
ments  of  these  lines  feed  into  a  range  computation  section  of  the 
tracking  lme  computer, 

8,  33  Comparison  of  Systems  One  and  Two 

It   is  obvious  from  comparing  Figs.    8.2,    8,3,    8.4  and  8.5, 
that  the  second  method  discussed  in  8.32  is  much  simpler  to 
design  and  construct  since  only  three  stellar  tracking  heads  are 
needed  while  the  system  discussed  in  8.  31  requires  six  tracking 
heads.     (Possibly  this  could  be  reduced  to  five  if  only  two  heads 
v/ere  used  for  space  reference  determination).     The  main  draw 
back  to  system  two  is  the  dependence  upon  gyros  for  space  refer- 
ence.    These  units  are  not,   at  present,    stable  enough  nor  pre- 
dictable to  the  degree  required     for  a  time  of  flight  of  the  length 
required  for  an  Earth -Mars  transit.     Section  8,64  discusses  this 
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SOLAR  SYSTEM 
TRACKERS 


ASSISTING  SPACE  REFERENCE 
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Fig.   8.4     Illustrative  schematic  diagram  of  tracking  system 
with  stellar  assisted  gyro  space  reference 
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particular  aspect , 

A  second  drawback  to  the  simplified  system  of  section  8„  32 
is  the  use  of  only  one  planet      The  effect  of  this  is  felt  in  com- 
putation of  (Rq.)  ,     j    and  will  be  treated  infection  8.41, 
8 .  4  The  Computing  Section 
8  41  Tracking  Line  Computation 

The  tracking  line  computation  section  of  the  Guidance  System 
receives  the  angular  measurements  of  the  tracking    lines  de- 
termined in  the  tracking  section.     Using  these  measurements  and 
data  from  the  Ephimeris  Storage  and  Time  Generator  within  the 
computing  section,    CR,0A^  u     anc*  (Vqa)  ,      are  determined „     In 
addition  the  various  separate  elements  of  the  ellipse  are  com- 
puted for  complete  identification  of  the  actual  trajectory  being 
followed  by  the  vehicle . 

The  method  of  computation  is  geometric  using  discrete 

measurements,    and  not  the  integration  of  velocities  and  acceler- 

71   72 
ations  experienced  by  the  vehicle,      '         The  reasons  for  this 

are  shown  in  section  8,6..     Chapter    7        developed  the  equations 

necessary  for  the  solution  of  this  problem  and  illustrated  the 

various  geometric  relations  in  Figs,  7,  2  through  7,  6, 

The  first  tracking  system  discussed  in  section  8,  31  has  as 
outputs  TLq  .    TL.0       and  TLq*  .     These  are  inputs  to  two  separate 
range  computers  which  determine  two  values  of  Rqa  using  the 
equation  of  section70  6  .     These  two  values  are  then  averaged  and 
smoothed  for  greater  accuracy  and  then  combined  with  the  vehicular 
heliocentric  latitude  and  longitude  (La,  Xa)  as  presented  by  TLq  , 
Since  a  series  of  Rqa  must  be  used  in  computation  of  the  elements 
of  the  elliptical  trajectory  ^   these  must  be  held  in  a  memory  storage 
unit,   and  fed  to  the  trajectory  elements  computer  upon  demand. 
The  determination  of  CVX.)  ,      is  made  by  a.  simple  integration 
process  and  along  with  (Rqa)   h     is  a  direct  output  of  the  tracking 
line  computer.     The  output  of  the  trajectory  elements  computer 
are    6,   a,  e^   and  i  observed.     The  simplified  information  flow 
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diagram  of  the  above  system  is  shown  in  Fig.    8,  6 

The  second  system  is  the  same  as  that  described  above  but 
uses  only  one  range  computer  with  the  resultant  loss  of  a  degree 
of  accuracy  in  the  solution,    although  a  system  simplification  is 
achieved  which  may  be  permissible.     A  simplified  information 
flow  diagram  of  this  is  illustrated  in  Fig     8   7 

8.42  The  Prediction  Computation  Section 


The  Prediction  Computation  Section  consists  primarily  of  a 
number  of  components  collectively  termed  the  Prediction  Com- 
puter,,   (see  Fig.    8   8).     The  purpose  of  the  Prediction  Computer 
is  to  continuously  establish,    as  a  function  of  time  after  t   .    the 
position  and  velocity  of  the  rocket  relative  to  the  reference  body 
The  position  and  velocity  so  established  may  be  termed  "pre- 
dicted" in  the  sense  that  they  do  not  stem  directly  from  obser- 
vations,  but  rather  are  values  which  should  actually  occur  if  no 
uncompensated  events  transpire  to  alter  the  rocket  motion.     The 
computer  must  establish  these  values  without  reference  to  any 
source  of  information  from  outside  the  rocket  itself.     This  in  no 
way  prohibits  the  computer  from  sharing  storage  data,    power 
supply,    time  generation,    and  the  like,    with  other  components. 
If  no  corrections  to  the  trajectory  were  envisaged  all  the  outputs 
could  be  computed  prior  to  launch     stored  in  a  memory,    and 
sampled  as  needed.     Since  corrections  are  envisaged  the  computer 
must  possess  the  capability  of  solving  values  for  the  rocket  po- 
sition and  velocity  during  the  entire  time  of  flight 

Selecting  the  Sun  as  the  center  of  our  inertial  reference 
frame  we  may  write  Equation  3-11 

ROA     j   =  fnf  +  G©A  +  (G®A  "  %q)  +  (G(T  A  '  G€o)  +  (G"dA  "  §&) 

(8-1) 
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From    Fig        2.8      we  determine 

G®©    <      10~'    ft/sec2 

GG©    <     10   9   ft/se°2       /=°  (8*2) 

Gcfo  <  10"8  ft/se°2 

Hence  the  equation  of  motion  to  be  solved  continuously  by  the 
Prediction  Computer  is 


R0A 


fnf  +  G©A  +  G©A  +  G<lA  +  GC?A  (8 '3) 


The  threshold  sensitivities  of  the  components  will  determine 
the  Phases  referred  to  in  Section  3    6.     Below  these  thresholds 
a  particular  component  will  be  inactive  and  the  problem  will  be 
reduced  to  the  appropriate  one  of  Equations   3-14,    3-15,   or  3-16. 

Referring  to  Fig     8   8  the  functions  of  the  sub -components 
are  almost  self -explanatory  when  examined  with  regard  to 
Equation  8  -3.     The  Summer  forms  the  quantity  representing  the 
acceleration  of  the  rocket  with  respect  to  the  Sun,     The  acceler- 
ation is  integrated  twice  to  form  the  rocket  position  vector, 
which  is  fed-back  to  obtain  values  for  those  gravity  fields  con- 
sidered in  the  equation.     Thrust  and  rocket  mass  form  the 
remaining  inputs  to  the  Summer 

8   43    Comparator  Section 

The  comparator  section  may  be  considered  an  integral  part 
of  the  main  computer      Its  function  is  to  compare  observed 
position  in  space,    (R(T)A)  -u    ,    emanating  from  the  tracking 
computer,    with  predicted  position  in  space,    (Rq   )pRF,.    derived 
from  the  prediction  computer.     Deviations  of  observed  position 
from  predicted  position  are  determined  and,  :  '  of  sufficient 
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magnitude,   form  the  basis  for  corrective  action.     The  test  of 
"sufficient  magnitude"  is  a  function  of  the  position  of  the  vehicle 
with  respect  to  Mars,     Obviously,   the  permissible  deviations 
would  be  small  during  the  early  stages  of  flight,   when  the  position 
of  the  rocket  relative  to  Mars  is  quite  lar^e.     Similarly,   as  this 
relative  position  becomes  less,   tolerances  become  larger  and 
errors  assume  a  lesser  importance,     A  deviation  computer 
weighs  or  assesses  the  rocket  position  relative  to  Mars  and 
creates  this  standard  of  measurement.     Fig.   8.9  is  a  functional 
representation  of  this  section. 

If  actual  deviations  in  position,    (D.)  R©a>    exceec*  those 
permissible  deviations,    (D  )  Rq,,   the  correction  comparator 
is  assumed  to  emit  a  signal  actuating  the  corrector  section. 
This  will  automatically  cause  initiation  of  the  computation  of 
the  new  trajectory. 

8.  44    The  Corrector  Section 

The  corrector  section  can  properly  be  considered  a  most 
complex  and  critical  part  of  the  computer.    Upon  its  proper 
performance  will  depend  an  accurate  approach  to  the  objective 
planet  or  the  alternative,   a  body  in  free  orbit  in  space  forever. 

The  corrector  section  is  continually  supplied  with  observed 
rocket  position,    (R-0A)  h   ,   and  the  elements  of  the  observed 
trajectory,   as  well  as  the  position  of  Mars  with  respect  to  the 
Sun  as  afunction  of  time,  R-Q^t).     These  quantities  form  the 
basis  for  the  computation  of  a  new  path  to  the  sought-after  planet. 
A  discussion  of  this  general  problem  has  appeared  in  Chapter  7. 

The  initiation  of  operations  within  this  section  will  be  de- 
pendent upon  receipt  of  the  comparator  section  actuating  signal, 
Sg(Act), 

The  modified  trajectory  is  assumed  chosen  through  an 
optimization  process  which  minimizes  fuel  expenditure,   time, 
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or  some  predesigned  parameter.     Once  this  selection  process 
is  completed     signals  proportional  to  the  required  velocity 
changes  can  be  converted  into  control  system  commands,    (S   )  T\ 
which  will  eventually  translate  the  changes  into  corrective 
action,  (See  Fig     8,10.  ) 

From  a  practical  viewpoint,    however,   the  selection  of  an 
alternate  path  in  flight  can  be  expected  to  pose  a  difficult 
problem      Size  and  weight  considerations  will  necessitate  an 
absolute  economy  in  functional  employment  of  components.     Even 
then,   the  optimization  problem  is  complex  and   may  well  repre- 
sent a  major  obstacle, 

8.  5  Thrust  and  Stabilization 

The  thrust  and  stabilization  unit  is  essentially  a  modification 
of  the  familiar  control  loop,   altered  because  of  operating    en- 
vironment to  include  rocket  thrust  device  rather  than  control 
surfaces,     A  signal  proportional  to  the  required  velocity  actuates 
the  thrust  device  which  is  positioned  prior  to  firing.     Several 
means  of  directional  control  are  available  for  corrective  thrust 
requirements : 

(1)  Movable  vanes  in  the  jet  stream 

(2)  Small,    individual  variable  thrust  jets 

(3)  Universally  mounted  main  jet  engine 

(4)  Gimbal  mounted  rocket  thrust  chamber 

Positioning  of  the  thrust  device  would  be  in  response  to  the 
correct  or  section  output  signal,   as  would  the  controlled 
metering  of  propel lant  to  insure  a  velocity  change  of  proper 
magnitude, 

The  accelerated  motion  of  the  rocket  produced  by  the  thrust 
device  would  then  be  sensed  by  an  accelerometer  which  com- 
pares    in  effect,    commanded  thrust,    with  resulting  thrust. 

The  stabilization  of  the  rocket  in  space  is  another  problem 


of  no  little  difficulty.     Some  of  the  proposals  for  rocket  orien- 
tation have  included  micro-jet  control,   gyroscopic  torque 
stabilization,   and  flywheels .     Many  authors  combine  vehicle 
motion  and  the  stabilization  although  they  distinctly  differ  and 
require  functional  separation. 

Fig,    8,11  is  a  line  functional  information  flow  diagram  of  a 
proposed  thrust  and  stabilization  sections. 

8.6  Errors  and  Effects 
— — — -r- 

It  has  been  stated  previously  that  the  unpowered  elliptical 
orbit  trajectory  for  transfer  between  planets  is  subject  to  large 
errors  in  terminal  position  caused  by  very  small  perturbing 
forces  during  transit  and  also  very  slight  errors  in  velocity 
vector  orientation  at  the  end  of  the  initial  thrust  period.     If,   as 
has  been  assumed  in  this  paper,    corrective  action  must  be  taken 
during  mid-course,   the  reduction  in  terminal  error  will  be  a 
function  of  time  remaining  till  intercept,    accuracy  in  determining 
the  actual  trajectory  being  traveled  and  the  accuracy  of  determi- 
nation and  orientation  of  the  corrective  velocity  vector.     We  will 
now  show  the  magnitude  of  some  of  the  terminal  errors  caused 
by  various  assumptions  and  inaccuracies . 

8.  61  The  Assumption  of  Circular  Coplanar  Orbits 

The  assumption  of  coplanar  circular  orbits  has  been  made 

by  numerous  authors  for  the  sake  of  simplicity  in  presenting 

.  .,  22.30,31,40,62,68        ,       .  „  ., 

basic  ideas.  and  others.    If  we  compare  the 

orbits  of  Earth  and  Mars  and  transfer  between  these  orbits  we 

find  the  following  to  be  true  for  the  worst  possible  cases. 

At  "favorable  opposition"  (once  every  16  years)  Mars  is 

fi 
35  x  10     miles  from  Earth.     At  "unfavorable  opposition"  Mars  is 

fi 
63  x  10     miles  from  Earth.     This  gives  a  maximum  actual  vari- 

ation  in  distance  between  orbits  of  28  x  10     miles  (.3  AU),   while 

the  so  called  average  distance  is  only  49  x  10     miles  (.  523  AU). 
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This  is  a  57%  variation  which  is  ignored  with  the  assumption  of 
circular  orbits.     Using  the  coplanar  assumption  (not  taking  the 
actual  orbital  inclination  of  1     51"  into  account)  we  introduce  a 
terminal  vectical  error  of  4.  57  x  10     miles  for  the  worst  possible 
case. 


8.  62 Errors  in  Astronomical  Distance 

The  basic  unit  of  astronomical  distance  is  the  Astonomical 
Unit  (AU).     This  is  determined  by  the  mean  range  of  the  sun  from 

the  Earth  and  by  international  agreement  (Paris,   1911)  has  been 

fi 
established  to  be  93  x  10     miles.     However,   this  distance  is 

based  on  angular  measurements  of  the  diameter  of  the  sun.     Due 

to  limits  of  resolution,   limited  base  lengths  for  triangulation 

and  other  effects,   this  basic  unit  is  thought  to  have  a  probable 

54   69 
error  of       10,  000  miles.      '  This  uncertainty  is  therefore 

carried  over  to  the  measurement  of  positions  of  solar  bodies 

and  hence  to  vehicular  orientation  in  solar  space.     Until  more 

accurate  astronomical  measurements  are  available  a  rather  large 

residual  error  in  spatial  position  will  have  to  be  accepted  and 

worked  with. 


8.  63Errors  in  Initial  Velocity  Vector 

If  we  now  take  a  typical  trajectory  from  Earth  to  Mars, 

we  can  see  the  effects  of  errors  in  initial  velocity  vector  on 

fi 
terminal  positions.     Assume  a  trajectory  of  arc  length  354  x  10 

miles  such  as  used  in  section  4.14   .     This  trajectory  has  a  time 
of  flight  of  172  days.,    and  an  orbital  inclination  of  2 .  51    .     An  initial 
velocity  error  of  1  foot  per  second  for  this  trajectory  will  give 
an  error  in  arc  length  alone  of  2  820  miles.     This  is  neglecting 
the  effects  on  the  actual  elliptical  properties  of  the  orbit  which 
increase  the  actual  miss  distance  even  more.     Using  an  average 
orbital  velocity  of  23.  8  miles  per  second,   this  would  be  equiva- 
lent to  an  error  in  arrival  of  118.5  seconds,   at  which  time  Mars 
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(traveling  approximately  15  miles  per  second)  would  be  1780  miles 
behind  the  initially  proposed  intercept  point. 

An  error  in  inclination  or  angular  orientation  of  the  initial 

velocity  vector  of  .  01  degrees  gives  a  terminal  error  normal  to 

the  trajectory  of  18,  600  miles.     We  can  also  approach  initial 

angular  error  from  the  standpoint  of  what  will  be  the  requirement 

for  insuring  a  terminal  error  of  not  more  than  1,  000  miles.     This 

will  necessitate  an  accuracy  within  1.  94  second  of  arc. 

49 
J.  G.   Porter       discusses  the  effects  of  small  initial  velocity 

vector  errors  on  the  shape  of  a  simple  Earth-Mars  elliptical 

trajectory.     This  trajectory  starts  at  its  perihelion  and  arrives 

at  Mars  at  aphelion  and  has  a  basic  time  of  flight  of  281  days.   Porter 

concludes  that  in  order  for  this  trajectory  to  arrive  at  Mars 

within  a  distance  of  50,  000  miles,  the  velocity  must  be  controlled 

to     001  miles  per  second  or  one  part  in  18,  000.     This  would  mean 

that  to  get  within  1,  000  miles  of  the  desired  terminal  position  the 

initial  error  would  have  to  be  less  than  .106  feet  per  second. 

For  this  same  trajectory  Porter  concludes  that  an  error 
of  12  seconds  of  arc  in  initial  velocity  vector  orientation  will  lead 
to  a  50,  000  mile  error  at  termination.     This  means  that  for  a 
1.  000  mile  error  we  would  have  to  be  within  ,  4  of  a  second  of 
arc  at  burn  out. 

The  major  conclusion  which  can  be  drawn  from  a  look  at 
errors  in  initial  velocity  vector  orientation  is  that  a  simple 
unpowered  trajectory  with  no  corrective  thrust  applied  during 
mid-course  will  not  suffice  for  interplanetary  travel.     This  was 
previously  stated  in  section, 4.  4      HDwever     of  the  two  errors  con- 
sidered,  velocity  and  angular,   the  errors  in  angular  orientation 
are  the  more  critical, 

8,  64  Mid  -Course  Errors 

We  can  now  look  at  errors  in  some  instruments  which  might 
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be  used  within  the  vehicle  for  guidance  purposes  such  as  star 
trackers,   accelerometers,   gyro  packages,   and  associated 
equipment. 

If  we  presume  the  need  for  a  maximum  star  tracker  error 
equal  to  or  less  than    1,  000  miles  in  determining  the  positions  of 
solar  bodies  1  AU  distant,   we  obtain  a  limit  of  2.22  seconds  or 
,0382  minutes  of  arc.     This  is  a  maximum  allowable  closed  loop 
error  composed  of  errors  in  reference  plane,   tracking  head, 
and  all  associated  equipment. 

Another  approach  to  this  problem  of  star  tracker  error  is 
to  assume  that  each  sighting  is  a  random  event  with  Gaussian 
distribution  of  the  error  in  angular  measurement.     If  we  now 
postulate  that  we  have  a  "maximum"  error  of  -5  seconds  of  arc 
in  the  values  obtained  from  a  star  tracker  and  that  99.  7%  of  the 
measurements  will  lie  within  these  values,    95%  will  be  in  error 
less  than  3,  33  seconds  while  68%  will  be  in  error  less  than    1.67 
seconds,     Since  an  error  of  1  second  of  arc  is  equivalent  to  450 
miles  at  a  distance  of  1  AU,   we  can  assume  that  68%  of  the 
measurements  would  have  an  error  less  than  750  miles,   and  95% 
would  have  an  error  less  than  1,  500  miles  when  sighting  bodies 
1  AU  distant.     It  appears,   then,     if    a   number  of  readings  are 
taken  over  a  short  period  of  time  with  proper  smoothing  and 
averaging,    suitably  accurate  values  can  be  obtained  from  tracker 
systems  which  may  by  themselves  have  a  considerable  possible 
error, 

To  determine  an  upper  limit  of  the  threshold  of  sensitivity 
to  acceleration  required  for  the  trajectory  of  Section  4.14  of  this 
paper  we  can  again  use  a  maximum  error  of  1,  000  miles  at 
termination.     This  corresponds  to  4.75  x  10       feet  per  second 

for  a  time  of  flight  of  172  days ,     For  this  same  trajectory  the 

71 
maximum  allowable  drift  of  a  purely  "inertial"      guidance  system 

to  insure  vertical  displacement  error  of  not  more  than  1.  000  miles 

-5 
would  be  1,64  x  10       mr  per  day.     This  is  determined  using  the 
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trajectory  arc  length  of  354  x  10     miles  and  time  of  flight  of  172 
days . 

The  primary  conclusion  that  can  be  drawn  from  the  above 
simplified  mid-course  guidance  system  errors  is  that  although 
completely  self  contained  inertial  guidance  systems  are  highly 
desirable,   the  required  accuracies  are  rather  high;  particularly 
if  return  trips  are  contemplated  without  rezeroing  the  system. 
However,    with  the  addition  of  star  tracking  units,    suitably 
accurate  guidance  can  be  obtained  which  is  not  limited  by  or 
dependent  upon  time  of  flight. 
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CHAPTER  9 

CONCLUSIONS  AND  RECOMMENDATIONS 
FOR  FURTHER  STUDY 

9.-1        Introduction 

The  contents  of  this  paper  lead  to  several  conclusions  and 
recommendations  of  interest  to  those  pursuing  the  problem  of 
interplanetary  travel,     Many  areas  included  in  this  thesis  are 
presented  with  a  view  toward  future  use  by  students  of  the  field 
and  lead  to  a  better  definition  of  the  problem.     The  bibliography 
represents  a  comprehensive  literature  search  of  available  un- 
classified related  material. 

9,  2       Conclusions 

(1)  During  the  time  of  flight  of  a  rocket  from  Earth  to 
Mars  the  center  of  the  Sun  may  be  considered  an 
inertial  point, 

(2)  During  a  major  portion  of  the  trajectory  the  problem 
may  be  considered  one  of  two  bodies,   with  some 

perturbations. 

(3)  Threshold  sensitivities  and  drift  uncertainties  of 
current  gyroscopic  mechanisms  preclude  the  use  of 
a  pure  inertial  navigation  system     at  this  time,   par  = 
ticularly  when  a  return  trip  is  contemplated „ 

(4)  During  the  time  of  flight  of  the  rocket  the  motion  of  the 
"fixed"  stars  is  so  small  that  they  provide  a  means  of 
establishing  a  space  reference  frame. 

(5)  Determination  of  rocket  present  position  m  space  within 
the  limits  proposed  in  this  paper     lies     within  the   capa 
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bilities  of  the  current  "state  of  the  art". 

(6)  Within  the  limits  of  the  assumptions  made,   the 
"combination"  trajectory^  fi  ,  an  initial  thrust  period 
followed  by  free  fall  flight;  modified,   if  necessary,  by 
corrective  thrust)  appears  to  offer  the  most  immedi- 
ately feasible  method  of  travel  to  Mars . 

(7)  Rocket  velocity  vector  orientation  appears  to  be  more 
critical  than  magnitude. 

22 

(8)  The  Hohmann  ellipse   '"  provides  the  minimum  required 

velocities  to  get  to  Mars  but  if  rocket  capabilities 

include  V.        (as  defined  within  this  paper)  of  10,  000  ft /sec 
inc  r  v      *  *  i 

or  more  a  wide  range  of  trajectories  is  available. 

(9)  The  Rocket  Plane  Concept  provides  a  means  of  reducing 
the  trajectory  to  a  two-dimensional  problem,    thereby 
greatly  simplifying  rocket  construction  requirements. 

It  is  emphasized  that  certain  assumptions  often  made 
in  the  literature  to  reduce  the  problem  to  two -dim  ens  ions 
are  unrealistic  when  viewed  in  the  light  of  actual 
condition.     The  Rocket  plane  concept  does  not  assume 
two -dimensions  but  selects  the  reference  plane  such 
that  deviations  from  this  reference  plane  are  negligibly 
small. 

(10)  On  a  kinetic  energy  basis  a  certain  optimum  departure 
orbit  exists  for  any  given  trajectory  to  Mars ,    Similarly 
an  optimum  arrival  orbit  exists  about  Mars .    Additionally 
these  optimum  orbit    radii  are  characterized  by  the 
property  that  for  a  relatively  small  increase  in  available 
velocity  a  satisfactory  departure  orbit  radius  may  be 
used  which  is  much  nearer  to  the  Earth  surface  than 

the  optimum  orbit  radius. 

(11)  Instrumentation  for  travel  between  planets  appears  to 

lie  within  the  capabilities  of  the  present  "state  of  the  art". 


3      Recommendations  for  Further  Study 

(1)  Digital  computer  trajectory  analysis  should  be 
instituted  covering: 

(a)  Required  parameters  for  all  trajectories  whose 
initial  rocket  velocity  lies  within  reasonable  limits. 

(b)  The  effect  of  variations  in  initial  conditions  on 
linear  miss  distance  at  Mars. 

(c)  Perturbative  effects  of  solar  system  masses  which 
were  ignored  in  this  paper  due  to  their  small 
magnitude. 

(2)  Additional  study  should  be  instituted  on  the  design  and 
construction  of  free-fall  space  reference  stabilization 
systems  suitable  for  interplanetary  travel. 

(3)  Studies  should  be  commenced  which  would  lead  to  the 
establishment  of  optimization  standards  for  trajectories 
to  Mars  on  the  basis  of  both  fuel  expenditure  and 
kinetic  energy. 

(4)  A  study  of  permissible  deviations  from  prescribed  path 
during  flight  allowing  terminal  accuracy  within  pre- 
scribed limits  should  be  instituted. 
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SYMBOLS 
DEFINITION  SUMMARY  1 


Symbol  Definition 

0  Sun 

C  Moon 

0  Earth 

Cf  Mars 

%  Jupiter 

Saturn 
Uranus 

^  Neptune 

E  Pluto 

O  Mercury 

Q  Venus 

<r  star 

a  Semi  major  axis  of  any  ellipse 

b  Semi  minor  axis  of  any  ellipse 

c  The  magnitude  of  the  effective 

exhaust  velocity  vector  measured 
with  respect  to  the  walls  of  the 
rocket  nozzle 

e  Eccentricity  of  any  ellipse: 

defined  by  the  relation 

b2 

e  =  1  -  — 2 ;  a  >  b 

a 
ff  Unit  field  forces  -  forces  per 

unit  mass  whose  existence  is  due 

to  the  presence  of  a  field 
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Symbol  Definition 

f  -  Unit  non  field  force  -  forces  per 

unit  mass  whose  existence  is  due 
to  all  causes  other  than  fields 

i  Inclination  of  orbit  -  the  angle 

measured  from  reference  plane 
to  plane  of  orbit,    in  a  plane 
perpendicular  to  the  line  of  nodes; 
generally  positive  where  measured 
in  a  direction  "above"  the  refer- 
ence plane 

j  A  variable  indicator  which  may 

take  on  identities  as  propounded 
within  the  text 

m.  Instantaneous  mass  of  rocket 

p  An  arbitrary  planet 

r.  Magnitude  of  radius  vector  from 

the  reference  focus  of  an  ellipse 
to  the  indicated  point  on  that 
ellipse,     Subscript  j  refers  to 
indicated  quantity 

r.  Radius  of  the  circular  orbit  about 

the  indicated  body 

t  The  independent  variable  time 

t.  Burning  time  -  that  period  of 

time  wherein  a  given  amount  of 
propellant  is  consumed 

t  The  time  at  which  an  imaginarv 

c  B  " 

rocket  would  have  arrived  at  Mars 
future  position  if  it  had  departed 

aphelion  of  its  ellipse  at  t  =  t 

Note:  t »  =  t      -  t . 
fed 
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Symbol  Definition 

t.  The  time  at  which  an  imaginary 

rocket  would  have  arrived  at 

Earth  present  position  if  it  had 

departed  aphelion  of  its  ellipse 

at  t  =  t 
o 

t-  Time  of  flight 


uf 


t.  Time  at     j     instant;  subscript 

j     refers  to  indicated  quantity 

t  The  instant  of  time  at  which  a 

o 

given  operation  commences. 

A  Angle 

AU  Astronomical  Unit  -  the  distance 

associated  with  the  mean  parallax 
of  the  sun.     By  international  a- 
greement  (Paris  1911),   this  was 
set  as  8.  80"  of  arc  and  corre- 
sponds to  about  93,  000,  000  miles, 
In  this  text  the  numerical  value  of 
the  AU  will  be  taken  as  93,  000,  000 
miles . 

D  Drag  -  as  conventionally  defined 

E  The  ecliptic  plane  -  the  plane 

defined  by  the  loci  of  points 
representing  the  instantaneous 
position  of  the  earth  throughout 
the  year  measured  with  respect 
to  the  sun 

£  Eccentric  Anomaly 

F  Force  -  as  conventionally  defined 
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Symbol  Definition 

G  Gaussian  gravitational  constant 

approximately  equal  to: 

-8       3  2 

6.67  x  10 "    cm    /gram  (m)  sec 

G.  Directed  acceleration  due  to 

J 

gravity  field  of  the  indicated 

body.     A  second  subscript  refers 

to  that  body  experiencing  this 

directed  acceleration.     Subscript 

"j"  refers  to  indicated  body 

GM.  The  product  of  the  Gaussian 

J 

gravitational  constant  and  the 
mass  of  the  indicated  body.   Sub- 
script "j"  refers  to  indicated 
body 

Imp  Impulse  -  the  integral    of  the 

magnitude  of  thrust  over  the 
time  of  its  duration 

KE  Kinetic  Energy  -  as  conventionally 

defined 

L.  Celestial  heliocentric  latitude  - 

3 

angle  measured  in  a  plane  perpen- 
dicular to  the  ecliptic  plane  from 
the  ecliptic  plane  to  the  line  drawn 
from  the  sun  to  the  indicated  point; 
measured  positively  toward   ce- 
lestial north 

L  Lift  -  as  conventionally  defined 

M  Total  mass  of  propellant 

M  Gross  mass  of  rocket  at  instant 

of  firing 
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NN'  Line  of  nodes  -  line  drawn  through 

the  ascending  and  descending  nodes 

PE  Potential  Energy  -  as  conventionally 

defined 

PR  Payload  Ratio  -  PR  =  1  -  £ 

R  Vector  from  the  reference  point 

to  the  indicated  point 

R  Magnitude  of  R 

T  Thrust  -  reaction  force  on  the 

rocket  or  vehicle  unit  structure 

TE  Total  Energy  -  as  conventionally 

defined 

V/tn  .  Velocity  vector  of  the  indicated 

body  with  respect  to  a  heliocentric 
non  rotating  coordinate  system. 

V.  .  Velocity  vector  of  the  rocket 

relative  to  the  indicated  body, 
and  measured  in  a  heliocentric 
non  rotating  coordinate  system 

V.  Velocity  vector  of  the  rocket  in 

circular  orbit  about  the  indicated 
body  measured  in  a  body  centered 
non  rotating  coordinate  system 

V.  Velocity  vector  of  "escape"  from 

the  gravity  field  of  the  indicated 
body,    equal  to 

\[T    x    V". 

js 

measured  in  a  body  centered  non 
rotating  coordinate  system 
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V/.  Av  Velocity  vector  of  the  rocket  with 

respect  to  the  indicated  body 
taking  into  account  the  kinetic 
energy  required  to'overcome  the 
effect  of  the  gravity  field  of  the 
indicated  body,    measured  in 
heliocentric  non  rotating  coordi- 
nate system 

V.  That  increment  of  V,.  Av     which 

mc  (jA)g 

must  be  provided  by  the  rocket 


a  Angle  measured  from  the  x-axis 

to  the  thrust  vector 

/3  Angle  measured  from  the  y-axis 

to  the  thrust  vector 

7  Angle  measured  from  the  z-axis 

to  the  thrust  vector 

x,  y,   and  z  axes  represent  an 

orthogonal  coordinate    set.     The 

thrust  vector  has  its  origin  at 

the  origin  of  this  coordinate  set 

M 
£  Propellant  loading  factor  =  M" 

n.  0-io    for  the  orbit  of  indicated 

body 

0  The  angle,   measured  in  the  rocket 

plane  of  motion,   between  a  body- 
sun  line  and  rocket  ellipse  aphelion- 
sun  line;  positive  angles  when  di- 
rection of  measurement  is  opposite 
to  rotation  of  planets  about  the  sun. 
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e 

at  any  time  j 

e  at .  t  =  t  , 
d 

e 

at  t  =  t 
c 

0 

of  Mars  at  t  = 

t 

Symbol  Definition 

6  (contd)  Prime  (')  indicates  measurement 

from  perihelion  vice  aphelion  with 
direction  of  measurement  positive 
when  same  as  rotation  of  planets 
about  sun 

6. 

J 

*d 

e 

c 

es  .  __  «„-  -d 

X.  *  Celestial  heliocentric  longitude 

-  angle  measured  in  the  ecliptic 
plane  from  the  reference  direction 
Aries  (  HP  h    in  the  direction  of 
rotation  of  planets  about  the  sun, 
to  the  projection  on  the  ecliptic 
plane  of  the  line  drawn  from  the 
sun  to  the  indicated  point « 

it .  Longitude  of  perihelion  of  in- 

dicated orbit 

</>  The  angle,   measured  in  the 

rocket  plane  of  motion  included 
between  the  velocity  vector  of 
the  rocket  or  vehicle  and  a  line 
perpendicular  to  the  radius  vector 
from  the  sun  to  the  vehicle 

0d  ^  at    t  =  td 

</>  <£  at    t  =  t 

ikqs A  The  angle  measured  in  the  rocket 

plane  of  motion  included  between 
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Symbol  Definition 

iL>.  (contd)  V©a  anc*  ^e  ProJecti°n  °f  ^Offi 

on  the  rocket  plane. 

4f~f a  The  angle  measured  in  the  rocket 

plane  of  motion,    included  between 
V/^.  and  the  projection  of  Vq^ 
on  the  rocket  plane. 

co.  The  angle  included  between 

perihelion  and  the  line  of  nodes 
measured  in  direction  of  ro- 
tation of  planets  for  indicated 
orbit. 

A  Rocket  or  vehicle 

<"p  The  period  -  elapsed  time  between 

repetitive  events 

The  vernal  Equinox  (the  first  point 
of  Aries)  defined  as  the  direction 
indicated  when  projection  of  the 
line  of  intersection  of  the  earth 
equatorial  plane  and  the  earth 
ecliptic  plane  passes  through 
the  center  of  the  sun;  this  action 
occuring  as  the  sun  apparently 
passes  from  south  to  north  of  the 
equator 

£2.  Longitude  of  ascending  node  of 

J 

indicated  plane 
y.  Longitude  of  descending  node  of 

J 

indicated  plane 

Ascending  node  -  celestial  heli- 
ocentric longitude  occupied  by  a 


Symbol  Definition 


body  passing  from  "below"  to 
"above"  the  ecliptic  plane. 
Descending  node  -  Ascending 
node  plus  180     "Above"  ecliptic 
defined  as  that  part  of  space 
wherein  an  observer  would  note 
counter-clockwise  rotation  of 
planet. 
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Fig.  4.12     Earth  angular  position  measured  from  Earth  aphelion 
as  a  function  of  time    (Overlay  #1) 
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